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Chapter 1 
Introduction 



These are the notes of a lecture given in the spring semester of 2008, entitled "optical prop- 
erties of semiconductors" . The goal of this lecture was to give a broad overview of the basic 
physical processes that govern the interaction between the light and semiconductor. My goal 
was to show the richness of the topic and to show the thread connecting the original research 
of the sixties and today's litterature. Whenever possible, I tried to show the connection be- 
tween the fundamental aspects and the applications in today's devices. At this stage, these 
lecture notes are rather patchy and are not trying to substitute for a complete text. 

1.1 bibliography 

1.1.1 Introductory texts 

Here I mention introductory textbooks that have a section on optical properties of semicon- 
ductors. 

C. Kittel "Introduction to solid state physics" . Chapts 1-3 treat elementary aspects 
of optical properties of solids 

Ashcroft & Mermin," Solid state physics" .(Saunders College) In this classic Solid State 
physics reference textbook, chapter 28 and 29 summarize the semiconductor general 
properties in an fairly elementary level. 

K. Seeger "Semiconductor Physics, an introduction" (Springer). Contains numer- 
ous derivations of analytical expressions. Chapt 11-13 treat optical properties. 

M. Balkanski and R.F. Wallis, "Semiconductor Physics and Applications" (Oxford). 
Very good treatement, includes some fairly advanced treatements of the optical prop- 
erties (Chapt 10,11,14,17) 

1.1.2 Advanced treatements 

This list contains the references I used to prepare the lecture. 
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E. Rosencher and B. Winter "Optoelectronics" . Very good presentation, more ori- 
ented towards devices. 

H. Haug and S.W. Koch "Quantum Theory of the Optical and Electronic Properties of 15 

(World Scientific). Advanced treatements of the semiconductor optical properties, 
more oriented towards theory. 

Yu and Cardona "Fundamentals of semiconductor" (Springer) Excellent reference book, 
geared more towards the fundamental aspects of the optical properties. 

G. Bastard "Wave mechanics applied to semiconductor heterostructures" (Les edi- 
tions de physique). Very good description of the computation of electronic state in 
semiconductor heterostructures. 

P. Zory "Quantum well lasers" (Wiley?) Contains some good chapters on quantum 
well semiconductor lasers, especially one written by S. Corzine. 



1.2 Notes and acknowledgements 

Under this form, this script is meant for "internal use" since it does not always give proper 
credit to the figures, especially for the ones taken from the books cited above. The author 
would like to acknowledge Giacomo Scalari for his help in preparing this lecture, as well 
as the one from Profs L. Degiorgi, Prof. B. Devaud-Pledran, Prof. A. Fiore who nicely 
volonteered with some of their data. Tobias Gresch prepared the artwork of the cover. 



Chapter 2 

Introduction to Semiconductors 



In this chapter, we would like to rapidely review the key concepts on semiconductors. See 
Chapter 28 of Ashcroft and Mermin's book. 



2.1 Crystalline structure and symmetries 

A perfect crystal is invariant under the translational symmetry 

f^f+ua + vb + wc (2.1.1) 

where u,v,w are integers. 

The crystalline structure of GaAs is displayed in Fig. 2.1, and is of the ZincBlende type. It 
consists of two interpenetrated diamond lattices. 




Figure 2.1: Zinc-blende crystalline structure 
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Figure 2.2: Zinc-blende crystalline structure 
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(2.1.2) 



+ 4 



per cell. 

The structure does not posess any center of symmetry. The basis of this zinc-blende structure 
are thetrahedrical bonds Si-Si or Ga-As. Each atom is the origin of a double bond towards its 
nearest neighbors, as shown in Fig. 2.3. The absence of inversion symmetry is also apparent 
on these tetrahedric bonds. 




Figure 2.3: Terahedrical bond of the zinc-blende 



2.1.1 Wigner-Seitz cell 

A Wigner-Seitz cell is the region of the space that is closer to a specific point of the lattice 
than to any other point of the lattice. A primitive Wigner-Seitz cell is constructed by taking 
the perpendicular bisector planes of the translation vectors from the chosen centre to the 
nearest equivalent lattice sites. 



2.2. WAVEFUNCTIONS OF THE CRYSTAL, BLOCH THEOREM 
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2.1.2 Reciprocal lattice 

— * 

The set of all wave vectors K n that yield plane waves with the periodicity of a given Bravais 
lattice is known as its reciprocal lattice. The bases vectors of the reciprocal lattice are given 
by: 

2n 

A = - bxc (2.1.3) 

a-(bxc) 

and by rotation of indices. The reciprocal space exhibits translational periodicity with all 
equivalent reciprocal wavevectors spanned by the set of integers h,k,l 

G = hA + kB + lC (2.1.4) 

— * 

The set of wavevectors G is a basis into which the crytal potential may be expanded: 

V(r) = ^V^exp(iGf) (2.1.5) 

G 

2.2 Wavefunctions of the crystal, Bloch Theorem 

The Hamiltonian of a semiconductor crytal has the translation symmetry 

H(f+R) = H(f) (2.2.6) 
with R being a reciprocal lattice vector. 

The Bloch theorem states that the wavefunctions have two "good" quantum numbers, the 
band index n and a reciprocal vector k such that the wavefunctions of the crystal may be 
written as: 

i; nk (r) = ^ r u nk {v) (2.2.7) 

where u n k(r + R) = u n k(r) exhibits the periodicity of the crystal. 
It may also be written as: 

^ nfe (r + R) = e jk - r ^(r). (2.2.8) 

2.2.1 Proof of Bloch's theorem 

Ashcroft & Mermin, chapt 8 
Translation operator 

T R ^(r) =V(r + R). (2.2.9) 
If R belongs to the Bravais lattice, the operator will commute with the Hamiltionian 

T R fty(r) = H(r + R)^(r + R) 

= #(r)V(r + R) = HT^{v) 1 ' ' ' 

and therefore 

T R H = HT K . (2.2.11) 
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These two operators form a common set of commuting observable, we can therefore write 
the wavefunctions using eigenfunctions of both: 

Hip = EipT^ip = c(R)ip. (2.2.12) 

Considerons maintenant les proprietes de c(R) 

TrTb,' = T R+R / (2.2.13) 

et done 

c(R)c(R') = c(R + R'). (2.2.14) 

To simplify somewhat the notation, a, with i = 1,2,3 the basis vectors of the primitive 
Bravais lattice, and bi the basis vector of the reciprocal lattice. We can always define 
(because c(aj) is always normalized to unity: 

c ( a .) = e 2wix \ (2.2.15) 

If now R is a translation of the Bravais lattice, it can be written as 

R = n\a\ + n^ai + 1130.3, (2.2.16) 

then, using2.2.15 

c(R) = c(niai + n 2 a 2 + n 3 a 3 ) 

= c(a 1 )" 1 c(a 2 )" 2 c(a 3 )" s 

i=1..3 

Using the orthogonality relations between the direct and reciprocal basis vectors, written as 

a; • bj = 2tt%. (2.2.18) 

Considering the product between the vector R of the real space and the vector k of the 
reciprocal space, written in the basis of the bi: 

k= J2 k ^ ( 2 - 2 - 19 ) 

i=1..3 

and then 

k-R = 27r(A;ini + A;2n2 + A;3n3). (2.2.20) 
It was in fact the result we had obtained for c(R), that was written as: 
c(R) = exp(27ri(x 1 n 1 + x 2 n 2 + x 3 n 3 ) 

= exp(ik • R) 

setting naturally X{ = hi. The wavefunctions may well be written as: 

T R iP = iP{r + R) 

= c(R)ip = e ik ' R ^(r) 

that is one of the equivalent formulations of the Bloch theorem. 



(2.2.17) 



(2.2.21) 



(2.2.22) 



2.3. ELECTRONIC STATES: FORMATION OF BANDS AND GAPS 
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2.3 Electronic states: formation of bands and gaps 

Bragg reflection and the perturbed electron 

2.4 Band structure of III-V and group IV semiconduc- 
tors 

Look at the band structure of Si, GaAs, InP 

2.4.1 Group IV semiconductors (Si,Ge) 

The band structure of these semiconductors is very similar because: 

1. They do crystallize in the same crystallographic structure (diamond) 

2. they have similar electronic outer orbitals 

The structure of silicon is purely covalent. The last orbital of atomic silicon has the electronic 
configuration 3s 2 p 2 . There are therefore 4 electrons (2s et 2p) sharing an orbital that could 
contain 8 (2 for the s orbital, 6 for the p orbital). Silicon has therefore 4 valence bands. The 
band structure of silicon and germanium, two most important semiconductors formed using 
the column IV of the periodic table, is shown in Fig. 2.4. 




WAVE VECTOR 

Figure 2.4: Germanium (left), Silicon (center) and Gallium Arsenide (right) band structures. 

The valence band maximum is at k = and is degenerate with the heavy and light hole 
bands. A third important valence band is the "spin-split" called this way because it is split 
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by the spin-orbit interaction. Finally, most important in the band structure of Silicon and 
Germanium is the fact that the minimum of the conduction band does not coincide with the 
maximum of the valence band. The semiconductor is called "indirect" (See Fig. 8.1. 

The conduction band minimum in silicon is in the direction [010] and, as a result, also in 
the directions [OTO], [001], [001], [100], [TOO] for a total of six minima. 
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Figure 2.5: Minima of the conduction band of Si, Ge and GaAs 



In Germanium, in contrast, the conduction band minimum is in the directions corresponding 
to the cube's diagonal, and we have therefore 8 conduction band minima. 



2.4.2 III-V Semiconductors (GaAs, InP, ..) 



The band structure of III-V semiconductors is similar since the tetrahedral bonds have the 
same structure as the ones in Silicon or Germanium. In fact, the missing electron of the 
group III with the electron configuration As 2 4p (for example Gallium) is provided by the 
column V element (for example Arsenic) of configuration 4s 2 4p 3 and these bonds have a low 
ionicity. 

For a large number of III-V semiconductors, the bandgap is direct. 



2.5. EFFECTIVE MASS APPROXIMATION 
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Figure 2.6: Structure de bande du GaAs 
As an example, the computed band structure of GaAs is shown in Fig. 2.6. 

2.5 Effective mass approximation 
2.5.1 Band extrema, effective mass 

Even though the band structure of semiconductors is in general very complex, as apparent 
in Fig. 2.6, the fact that the Fermi level lies in the middle of the forbidden bandgap enables 
us to consider the band structure in a region close to the band extrema. In that case, the 
band dispersion can be expanded e(k) into a second order Taylor series: 

e(k) = e + V |^ • fa - k ) 2 . (2.5.23) 
i=i..3 aK * 

The equation 2.5.23 is written in the simplified case that the energy ellipsoid are aligned 
with the x, y and z axis of the coordinate system; in the general case terms of the form d ^.g k . 
must be added. 

By analogy with the free electron case, the term of the expansion can be identified with the 
inverse of an effective mass: 

J- - ww- < 2 - 5 - 24 ' 

'" / x,y,z 11 uh/ x,y,z 
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2.5.2 Semiclassical electron dynamics 

The analogy may be pushed further: the group velocity of the electron is: 

1 d(huj) lde 

and the electron will obey the semiclassical electron dynamics given by: 

F = h^- (2.5.26) 
at 

where 

F = —qS — qv x B (2.5.27) 

is the classical force on the electron 

2.5.3 The hole 

The motion of the electrons in a not completely filled band with negative mass can be 
interpreted with a particle called the hole, that has the following properties: 

• Its mass is positive and given by = — ^fp 

• Its wavevector is the opposite of the one of the missing electron: kh = —k e 

• Its charge is positive and the opposite of the one of the electron qh = —q e 

2.6 Density of states in 3, 2 and 1 dimensions 
2.6.1 3D 

Wavefunction 

ip(x, y, z) = Ae i{kxX+k v yk * z) . (2.6.28) 

In a finite volume, Born- von Karman periodic boundary conditions, ip(x+L, y, z) = if)(x, y, z) 
imply: 

k x = n x 2n/L 

k y = n y 2n/L (2.6.29) 
k z = n z 2n/L 

Using for the energy: 

h 2 k 2 

E = (2.6.30) 

and the number of state as a function of k m nr 



2.7. PHONONS 
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we have 



Differentiating with respect to the energy: 



2.6.2 2D 



2.6.3 ID 



^ /r ^ dN(E) 32 1 ,2m* r- 1 (2m*) 3 / 2 r- 



D(E) = (2.6.34) 



2.7 Phonons 

optical phonon, acoustic phonon 

2.8 Doping 

Impurities may act either as deep level, shallow levels or alloy. 
2.8.1 Hydrogenoid donors 

Same energy spectrum as the Hydrogen atom but renormalized by the dielectric constant 
and the effective mass 

2.9 Carrier statistics 
2.9.1 Ferm-Dirac Statistics 

The probability of finding an electron a energy E in a bath at temperature T is: 

f(E,f*,T)= Tg^r— (2.9.36) 

eX Pl kT ) + 1 

where k is Boltzmann's constant and \i is the chemical potential. 
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Figure 2.7: Fermi-Dirac distribution for various temperatures 



The total number of electrons in the structure is the integral of the product of the density 
of state and the Fermi-Dirac distribution function. 



tot 



OO j 

D(E) 7W —, 

E=0 ; eX p(^)+l 



dE. 



(2.9.37) 



2.9.2 Holes 

Distribution of holes: for the holes, one applies the rule (hole = no electrons) and then 



f h (E,v,kT) = l-f e (E,ti,kT) 



( E-n 
\ kT 



) 



(2.9.38) 



2.9.3 Classical limit 

In the case where E— /x >> kT, the Fermi distribution can be simplified because exp(-^^) >> 
1 and then 



f(E,»,kT)*e W (^-). 



(2.9.39) 



This corresponds to a probability of occupation much below unity and therefore is called the 
classical limit. 



2.9. CARRIER STATISTICS 
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Figure 2.8: Regimes in which the classical limit can be used, a) High energies in a metal, 
b) Low doped semiconductor 



In this approximation, the total number of electron can be computed in the bulk case: 

1 (2m*) 3 / 2 r= — , E-a, 



n 



n 



1 (2m*) 3 / 2 



2tt 2 



h 3 



kT 
E 



)dE 



exp(fikT) / oojE - -E G exp(-— )dE. 
Je g kl 



2tt 2 h? 

To carry out the integrate 2.9.40, we use a change of variable: 

E — Eq 



X = 777 , dx = T77.dE 

kT kT 



1 (2m*fcT) 3 / 2 , 
n = — exp ( 



2tt 2 fr> ' v kT 

The last integral is performed using Cauchy's formulai: 



A* — Eq\ I i / ■■) , , , 
J / x 1 exp(— x)dx 



x 1 / 2 exp(— x)dx 



x=0 



On obtains then finally: 



n 



2( 



m e k B T, 3/2 



(2.9.40) 



f 00 E f°° E 

y/E - E G exp(-—)dE = Vkf J ^ 2 exp(-x)exp{—^;)dxkT ( 2 . 9 . 41 ) 



(2.9.42) 
(2.9.43) 
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The same approach can be carried over for the holes: 

A«exp(-^^). (2.9.44) 
Assuming that the zero of energy is on the top of the valence band: 

P = J^D h (E)f h {E)dE = 2^^f 2 e Xp (^). (2.9.45) 

The product between electron and hole concentration yields: 

P = 4 (SW,) 3/2 exp(-:^). (2 . 9 .46) 



n 



The produt n ■ p depends only on the crystal nature and its temperature. 
The number of electrons in the band is also written in terms of the quantum concentration 
for the conduction n c and the valence band n v (also called effective density of states in some 
texts): 

n = n c exp( ) 



p = n v exp(-— ) 

where 



kT ' (2.9.47) 



f m e kT, 3 /2 
,m h kT,3/2 



(2.9.48) 



depend on temperature T and effective mass m*. In virtually all IIII-V semiconductors, the 
valence band consists of a heavy hole and a light hole band, the concentration for each band 
must be summed together: : 



V = J D hh f h + Di h fi h 



^2.9.49) 



that can be written as 2.9.45 if one sets: 

™T = m 3 h { 2 + mf (2.9.50) 

which is called "density of state effective mass). 
Note that the product n ■ p is constantt: 

np — n c n v exp^—Eo/kT^ — n { . (2.9.51) 

rii is called the intrinsic concentration and is: 

E 

m = v /n ^exp(-^ ; ). (2.9.52) 



2.9. CARRIER STATISTICS 
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It is the thermal concentration of carriers in undoped semiconductors. 




Figure 2.9: Intrinsic concentration for Ge, Si, et GaAs as function of reciprocal temperature 
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Chapter 3 

Light-matter interaction 



3.1 Oscillator model 

Solving the equation for a (classical) harmonic motion with damping: 

m x + 2m 7X + m u x = —q£(t) 
yields a polarization P{uj) at angular frequency uj given by: 



n q 2 1_ 

m uj 2 + 2i^uj — uj'q 



Remembering the definition of the susceptibility as: 

P(u) = e x(Lu)£H 

the latter is then 

n g 2 1 



m e c^ 2 + 227a; — Uq 
The above equations can be rewritten in terms of the sum of two poles 

,2 r 1 1 



m LV 
with c^o = a/c<Jq — 7 2 

3.2 Dielectric function 



— c<; + 27 uj — uj + 17 



The dielectric function can be then derived as: 

,2 

jj^uo - uo + r: jj + .7. + r ' 
17 



e (l-^r( \ \ )) 

V 2u; n uj — uj n -\- 17 U + U n + 17 ' 
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where the plasma frequency. 



The dielectric function can be written separately for the real and imaginary parts (neglecting 
the negative frequency pole): 



x' = 4 - 1 



UJ p 10 — uJq 

2u (u; - lu ) 2 + 7 2 



and for the imaginary part: 



// // 
X = e r = — 
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Auj (uj - lo ) 2 + 7 2 



(3.2.7) 



(3.2.8) 



Lorenzian model 
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Figure 3.1: Real and imaginary part of the susceptibility in the Lorenzian model. 

3.3 Effective medium approximation 

Taking into account the effect of other resonances in the susceptibility: 



e = too I 1 



e < , 1 



2u' uj — oj' + i'j uj + uj' + ij 



^)) 



(3.3.9) 



The only trick is to remember to add the polarization. Because of the form of Equ. 3.3.9, 
one is then lead to redefine the plasma frequency as: 



OJr. 



n q z 



(3.3.10) 



3.3. EFFECTIVE MEDIUM APPROXIMATION 
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Example of "classical" dynamics and plasmon oscillations: subpicosecond plasmon excita- 
tions in doped semiconductors (Kersting, C. Unterrainer et al., PRL 1997). As shown in 
figure 3.2, the excitation is done by a femtosecond Ti:sapphire laser. 




Figure 3.2: Ultrafast excitation of doped semiconductor. 



Electron-hole pairs are created at the surface, and excite the plasmons of the doped layer. 




Figure 3.3: Schematic drawing of the sample. 



As expected, oscillations are observed at the plasma frequency uu. 
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Figure 3.4: Spectrum of the oscillations as a function of doping. 



For the example above, using n = 1.710 16 cm 3 , = 13.1, an effective mass of 0.065, a 
computed frequency of 1.2THz. is obtained. 



3.4 Kramers-Konig relations 

The susceptibility nas some important mathemtical features because it is a linear re- 
sponse function. The polarization at time t can be expressed as a function of electric field 
at past times through: 

/t poo 
x (t-t')£(t')dt'= / x (r)£(t-T)dT (3.4.11) 
oo JO 

where, because of causality, xij) — Oi/r < xi 00 ) is therefore analytic in the upper half 
plane. Using a Cauchy integral argument, it follows that 

XH = P 1^ . , X{V)) Tv (3-4.12) 

The latter equation implies a relationship between real and imaginary parts of the suscepti- 
bility, called the Kramers- Kronig relations: 

xV) = p rwM (3 . 4 . 13) 

t {u) = J± P r*JpL ,3.4,4) 

7T In V 1 - UJ 2 



3.4. KRAMERS-KONIG RELATIONS 
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Figure 3.5: Comparison between the real part of the susceptilibity computed by a Kramers- 
Kronig transform of the imaginary part and by the expression given by 3.2.8. The difference 
arises from neglecting the negative frequency pole. 



3.4.1 Sellmeir's equation 

Approximating the imaginary part of the susceptibility by 



and using Kramers-Konig relations, we obtain for the real part of the susceptibility: 



which predicts the dispersion far away from the resonance. A simple generalization of this 
equation is referred as Sellmeir's equation for the refractive index, usually written as a 
function of the wavelength instead of the angular frequency: 




(3.4.15) 




(3.4.16) 




(3.4.17) 
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Q2 0.4 0.6 tl.8 L2 1.4 l.ft 



Figure 3.6: Comparison between a Sellmeir's equation with three resonances (two in the UV, 
one in the mid-infrared) and the experimental refractive index for BK7 glass. 

3.5 Interaction between light and a quantum system 

Depending on the result of the microscopic model, there are many ways to introduce the 
optical response. 

3.5.1 Fermi's golden rule: a loss term to a propagating wave 

Let us assume the interaction Hamiltonian 

Hint = —qr- S sm(ut) (3.5.18) 
the scattering rate, using Fermi's golden rule, writes: 

R = \ = ^ \^\ H ^t\f)p{E f -Ei- hw). (3.5.19) 

The absorption of the electromagnetic wave by the quantum system is responsible for a decay 
of the latter with an absorption coefficient a, where 

I(x) = I exp(-ax) (3.5.20) 

The energy loss by the wave per unit volume is 

P/v = la (3.5.21) 

and the intensity is related to the electric field by 

I = -tQn re f r c£ 2 (3.5.22) 



3.5. INTERACTION BETWEEN LIGHT AND A QUANTUM SYSTEM 
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where n re f r is the refractive index and c the light velocity. Balancing the energy loss of the 
EM field with the energy gained by the quantum system yields: 

a = — (3.5.23) 

e n re f r ct z 

therefore: 

a = _nu^_ J2\(i\D\f)fp(E f -Ei-tuv) (3.5.24) 

3.5.2 A polarization field 

— * — * 

As a polarization field P reacting to the incident field E. The polarization is computed by 
evaluating the dipole operator in the time-dependent wavefunction: 

P(t)=n Q ^(t)\qx\iP(t)) (3.5.25) 

and using a first order, time dependent perturbation expansion for the wavefunction ip(t), 
the susceptibility x(u;) is 

xM = -nr £ l*"'l 2 ( i 1 i • L ^-) ( 3 - 5 - 26 ) 

ft ^ X U + Ui m + VJ LU-UJim + l^' 

that can also be rewritten using the oscillator strength, defined as: 

fmi = -^-\x m i\ 2 uj m i (3.5.27) 
in a fashion very reminiscent of the classical expression: 

XH = -— T — ( ) (3.5.28) 

3.5.3 A current density j 

As a current density j that is added to Maxwell's equations, in a crystal that has already a 
dielectric response e(cu): 

- - 3D 

V x F = j + — =oE- iueE (3.5.29) 
The equivalent dielectric function is written as: 

e ;(^) = ^M + e^) (3.5.30) 
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3.6 Momentum p and dipole z matrix elements: sum 
rule 

3.6.1 Relation between p and z matrix elements 



As the commutator between p and z is [z,p] = ih, applying it to the kinetic term of the 
Hamiltonian one can derive the relation: 



(<j>n\p\<fim) = im uj nm {(f) n \z\(f) m ) (3.6.31) 



3.6.2 Sum rule 

Taking the completeness of the (eigen) states 



X>n)(0n|=l (3-6.32) 

n 



and the relationship between position and momentum matrix element, the sum rule for the 
oscillator strength can be derived: 

J]/ni = l (3.6.33) 

n 



The sum rule appears in a number of context. It can be expressed as the integral of the 
absorption strength. In the next graph, intersubband absorption for three different samples 
with various energy level ladder are compared. As expected, the integrated absorption is 
constant within experimental error. 



3.7. DIPOLE OR MOMENTUM MATRIX ELEMENTS 
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Figure 3.7: Comparison of the intersubband absorption for various structures, with energy- 
levels schematically drawn close to the curves 



3.7 Dipole or momentum matrix elements 

Hamiltonian derived from 

H = (P ~ qA)2 (3.7.34) 
2m 

Use Coulomb gauge (VA = 0). For low intensity (neglect the term in A 2 ), we obtain as an 
interaction Hamiltonian: 

Ei = — —A- P. (3.7.35) 
m 

Because of the large difference between the light wavelength and the atomic dimension, the 
spatial dependence of A(r) is neglected inside the matrix elements. This is called the dipole 
approximation. The ratio of the spatial frequency of the light wave k p h to the electronic 
part k e is 



kph _ / Eph 



k e y E rest 2(m*/m ) 



(3.7.36) 



where E rest = 550keV is the rest mass of the electron. This ratio is ~ 10 5 for leV photons. 
The form commonly used is then: 

(V*|ffi|V>i> = ^^mr)\Pmr)) (3-7.37) 
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For a plane wave, A is paralell to £, for a wave polarized along z and propagating in the y 
direction: 

A(f, t) = A e z e i{ - ky ~^ + A* e z e- i{ky - wt) . (3.7.38) 

It is convient to use A being pure imaginary such both electric and magnetic fields are real: 
that yield £ = 2icjA and B = 2ikA . 



Chapter 4 

Optical properties of semiconductors 



4.1 Reflectivity measurements, ellipsometry 
4.1.1 Reflectivity 

Use a Kramers-Kronig relation to express the phase of the reflectivity from the spectrum of 
the intensity reflectivity: 

ln{p{uj')duj' 



oh = r 

n Jo 



uo' 2 — to 2 



(4-1-12 



4.1.2 Ellipsometry 

There, the dielectric constant is obtained from a measurement of the ratio of the p to s 
reflectivity o = ^ through: 

e(w) = sin 2 + sin 2 0tan 2 0^— -^] 2 (4.1.2) 




Figure 4.1: Schematic diagramm of an ellipsometer. 
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Figure 4.2: Dielectric function of GaAs. Note the influence of the surface oxide. 



4.1. REFLECTIVITY MEASUREMENTS, ELLIPSOMETRY 
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Figure 4.3: Dielectric function of GaAs including the gap. 
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4.2 Critical points and band structure, Van Hove sin- 
gularities 
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Figure 4.4: Dielectric function of Ge. compared with the computation. 




Figure 4.5: Germanium bandstructure with the critical points 



4.2. CRITICAL POINTS AND BAND STRUCTURE, VAN HOVE SINGULARITIES 
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Let us use the Fermi golden rule, and assuming the invariance of the matrix element with 
energy and wavevector. The transition rate is then given by: 



R " ° 1 i \V k {Ec-E^\t Ev=hjJ (42 - 3) 



where dS is the element of surface in k space defined by the equation 



E c (k) - E v (k) = hw (4.2.4) 



and C is a prefactor given by 



Anhq 2 2 



C = —^Ai (4.2.5) 



m 



Expanding the energy around a critical point in a Taylor expansion: 



3 



E c - E v = E + a i(h - k 0i) 2 (4-2.6) 



i=i 



enables a classification of the Van Hove singularities depending on the relative signs of the 
coefficients ac M and M 3 if they have the same sign (positive, resp negative), Mi and M 2 
for saddle points. The shape of the absorption is shown in Fig. 4.6. 

As an example, for a three dimensional crystal with the Mo, J dS = 4nk 2 and \S7k(E c — 
E v )\Ec-E v =hw = — yields an absorption coefficient proportional to 



{m r f/ 2 ^E^W 



(4.2.7) 
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Figure 4.6: Van Hove singularities 



The square-root shape of the absorption at the band edge is compared to the experiment for 
InSb in Fig 4.7. 



4.3. REFRACTIVE INDEX 



33 



II) 4 



> 

J? 



lo- 



in 



InSb 




• Experiments 

ahv=A(h(0-E g ) U2 

Corrected for variation 

of matrix element 



JL 



0.2 0.4 0.6 

Photon energy hv [eV | 



0.8 



Figure 4.7: InSb absorption, compared with the computations assuming constant and non- 
constant matrix elements 



4.3 Refractive index 



For many devices, an accurate knowledge of the refractive index and its frequency depen- 
dence is very important. The Kramers Kronig relations, together with a modelisation of the 
imaginary part of the dielectric function, allow a phenomenological model of the refractive 
index dispersion of III-V direct semiconductor to be derived. Assuming the imaginary part 
of the dielectric function e" with the form: 



e"{E) = r)E A 



(4.3.8) 



between the bandgap energy Eq and a final cutoff energy Ef and zero elsewhere, as shown 
in Fig4.8 
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Figure 4.8: Imaginay part of the dielectric function for GaAs and the Aframowitz approxi- 
mation 



Because the dielectric response of the III-V materials cannot be easily represented by a finite 
sum of oscillators, this is a better approximation than a "normal" Sellmeir's equation. From 
a integration of the Kramers- Kronig equation, the refractive index is obtained as: 



n\E) = 1 + ±[E} - E%) + ^E\E) - El) + (4.3.9) 



where r], Ef and Eq fully characterize the material. In the litterature, the parametes are 
given in terms of parameters of a Sellmeir's equation 

e'(E) = l + ^^- 2 . (4.3.10) 
These parameters are related to the ones of our model through 



E s = pEl e -El (4.3.11) 

and 

^m^wr (43 ' 12) 

The experimental results for a few relevant III-V materials are shown in Fig. 4.9 and 



4.3. REFRACTIVE INDEX 
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Figure 4.9: Computed dispersion of the refractive index 
the relevant material parameters in Fig. 4.10. 
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Figure 4.10: Energy gaps for various III-V materials for the refrative index 



36 CHAPTER 4. OPTICAL PROPERTIES OF SEMICONDUCTORS 



Chapter 5 



Bulk semiconductors: bandstructure 
and fundamental gap 

5.1 k p approximation 
5.1.1 Basic approximations 

Behind the effective mass approximation lies a very powerful approach to the computation of 
the band structure. It relies on the knowledge of the band structure at k = and expanding 
the wavefunctions in this basis. The Schroedinger equation for a crystal writes: 

2 fc2 

(2^ + V{r) + 4mp ((? X ' = E ^ (5 - L1) 

For simplicity, we drop the spin-orbit coupling term. The latter arises as a relativistic term: 
the motion of the electon in the field of the ion, the latter sees an equivalent magnetic field 
that operates on the angular momentum variable. Let us first compute the action of p on 
if), written in terms of Bloch wavefunctions so that 



lUO = ^X*(^ (5-1-2) 



we obtain: 



p- {e^ujf)) = -ihV ■ (e ikr \^) (5.1.3) 

= hke^u^T) + e&p- u n% {r) (5.1.4) 

= ^{p + hk)u n% {f). (5.1.5) 
Using the above relation, the Schrodinger equation is obtained for u n ^(f): 

+ + "' (F)Ks(F) = E " kU - em (5X6) 

The Hamiltonian H — H + W(k) may be splitted into a k-independent 

H Q = -£- + Vlf) (5.1.7) 
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and k-dependent part 

ft 2 k 2 ft -> 

W(k) = —— + —k ■ p. (5.1.8) 
2m m 

The solution of the the equation 

H u n0 (r) = E n0 u n0 (r) (5.1.9) 

are the energies of the band structure at the T point k = 0. The fundamental idea of the 
k ■ p approximation is to use the u n o(f) as a basis for the expansion of the wavefunction and 
energies at finite k value. In the simplest cases, taking an interband transition across the 
gap and looking at the conduction band, taking the second-order perturbation expansion: 

In the lowest order approximation, all other bands except for the valence band may be 
neglected, in which case the dispersion can be written as (taking the zero energy at the top 
of the valence band): 

Defining the Kane energy Ep = 2moP 2 such that 

2 

E P = — |(w c ,obk,o)| 2 (5.1.12) 
m 

the dispersion of the conduction band can be written as: 

h 2 k 2 , E P 

= ^ c + — 1 + 7 
2m i 

We then obtain the effective mass as: 



E ' {k) = E - + ^ 1 + t } (5 - L13) 



(*.*)-' = (mo)- 1 + (5.1.14) 

The Kane energy is much larger than the gap Ep >> Eg and is rather constant across the 
III-V semiconductors. As a result, the effective mass in inversly proportional to the band 
gap. 

One should be careful that some authors use the definition 

E P = 2m P 2 (5.1.15) 
(Bastard, for example), while others (Rosencher) use 

E P = P 2 . (5.1.16) 
At least they all use Ep having the dimension of energy!. 



5.1. K P APPROXIMATION 
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Figure 5.1: Fundamental gap Eq, Split-off energy A, conduction band effective mass and 
Kane energy Ep for various III-V semiconductors 



In the above expressions, the spin-orbit term can be included by replacing the operator p by 

h_ 



n = p+ ir —^(axVV). (5.1.17) 



5.1.2 Beyond the perturbation expansion 

A very powerful procedure is to expand formally the solutions of Equ. 5.1.6 in the solutions 
at k = 0, writing formally: 

u nk (r) =^2(&(k)um l0 (r) (5.1.18) 

rn 

and restricting the sum to a limited relevant subset of bands. Improved accuracy can be 
achieved by introducing more bands. In this basis, and projecting the equation onto the 
state u Mfli the Hamilton equation is written as: 

{Emo + ^ - K(k))c%\k) + J2 H Mm<£\k) = (5.1.19) 

where the kp Hamiltonian is 

h 

H Mm = k ■ (UM,o\p\ U m,0) (5.1.20) 

m,Q 

To the extend the matrix elements are known, equation 5.1.19 can be solved. 

5.1.3 Example: a two-band Kane model 

Let us write explicitely a two-band Kane model. The latter can be a fairly realistic model if 
one is only interested at the effect of the valence band onto the conduction one, replacing the 
three spin-degenerate valence bands (heavy hole, light hole, split-off) by an effective valence 
band. The u n k can then be expressed as: 

Unk = a c u c0 + a v u v0 (5.1.21) 

Replacing this expansion into Eq. 5.1.19, we obtain the following matrix equation: 



E + pL J^k- p cv 

L 2mo mo 

-S-jfe • fiL E v + P^ 

mo ^ cv 2mo 









V a v , 




\ a v J 



(5.1.22) 



CHAPTER 5. BULK SEMICONDUCTORS: BANDSTRUCTURE AND FUNDAMENTAL 
40 GAP 



Assuming E c — 0, E v — —E G , the solution of the matrix equation satisfies: 



h 2 k 2 n 2 h 2 n 2 



This second order equation in E can be of course solved directly; it is however more instructive 

mo 1 



to write it under the form a pseudo effective mass equation (remember Ep = -^-p 2 cv '- 



2rriQ E + Eq K J 



For k — > the above expression reduces itself to the result of the perturbation expansion: 

(m*)- 1 = (m )- 1 (l + ^). (5.1.25) 
Equation 5.1.25 can be expressed in a somewhat simplified form: 



m*(E)=m*(0)(l + ^-) (5.1.26) 
Eg 



that is commonly used in the literature. 



5.1.4 Realistic model 

To have a realistic description of the band strucure, the following steps are taken. 

In a first step, the valence band k=0 Bloch basis function are build using orbitals \X, Y, Z ||), 
in such a way as to build states that are eigenstates of the total angular momentum operator 
J 2 and J z . Of course at this point the direction of quantification is totally arbitrary. The 
figure 5.3 gives this expansion explicitly. This expansion is done in two steps: first one 
should eigenstates of the orbital angular momentum as a function of the orbitals and then 
use the sum of angular momentum rules to get the total angular momentum (orbital plus 
spin). 
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Figure 5.2: Schematic band structure 
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Figure 5.3: Definition of the basis vector 



Once this basis is defined, the matrix elements are taken between these basis states dropping 
the k-dependence of the spin-orbit term. Relativistic term such as the Darwin term are also 
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neglected. The only non-zero matrix elements between orbitals are 

P = — (S\p x \X) = —(S\ Py \Y) = —(S\ Pz \Z) (5.1.27) 
m m m 

between orbitals with the same spin because all others vanish by symmetry. The resulting 
matrix is shown in Fig 5.4. 
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Figure 5.4: Matrix equation for the k ■ p 8 band model 

Solving the matrix equation shown in Fig. 5.4, one obtain an implicit equation for the 
dispersion given by: 

A(jfe) = -e G (5.1.28) 
2A 

X(k) [X(k) + e G ] [X(k) + e G + A} = h 2 k 2 P 2 [X(k) + e G + — } (5.1.29) 
where X(k) is defined by 

A « = £ «-|C (5L30) 
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Equation 5.1.29 is the one of the heavy particule (at this level, they display a positive (i.e. 
electron-like) dispersion, and the equation 5.1.29 the light particules. The computed bands, 
with GaAs band parameters, are shown in Fig. 5.5. 
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Figure 5.5: Computed band structure of the GaAs using the 8 band model 
The structure of the valence band is shown in more detail in Fig. 5.6. 
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Figure 5.6: Computed valence band structure of the GaAs using the 8 band model 

A striking feature is the electron-like dispersion, with a bare electron mass, of the heavy 
hole states. It means that the heavy hole band is not coupled to the conduction band. 
The correct dispersion of the heavy hole band can only be predicted when including the 



GaAs 
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effects of the remote bands. The absence of coupling is apparent from inspection of the k • p 
by remembering that the direction of quantification is arbitrary and therefore can be chose 
always parallel to the vector k. In that case, k = k z and it is apparent that the states ||, ±|) 
are not coupled to any other band. For an other direction of quantification, the same heavy 
particle is a linear combination of ||, ±|) and ||, ±|) states. 

Effective masses can be obtained from this model through a Taylor expansion of the solution 
of equation 5.1.29 around k = 0. The values obtained are 



1 4P 2 
m 3 E G 


2 P 2 
+ 3 E G + A 


(5.1.31) 


1 

™r 8 


1 4P 2 
m 3 E G 


(5.1.32) 


1 1 


2 P 2 


(5.1.33) 
(5.1.34) 


m Vl m 


3 Eg + A 



where the group theory notation has been used to distinguish the various band edges: r 6 
corresponds to the conduction band, r 8 the light holes and the split-off band. The mirror 
effect of the bands is clearly apparent. 

Comparison between this dispersion and the results of a pseudo-potential computation is 
shown in Fig. 5.7 




Figure 5.7: Comparison between the kp approximation and a pseudo-potential computation 

The k ■ p approximation is very good as long as the wavevector remains close to the center 
of the Brillouin zone. In particular, it is unable to predict the minima at the edge of the 
Brillouin zone. 
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5.2 Computation of the absorption edge in bulk mate- 
rials 

From the comparison between the energy lost by the optical wave and the one gained by the 
quantum system, we had derived in equation 3.5.23 a relationship between the scattering 
rate given by Fermi's golden rule R and the absorbtion for an electric field S. 

2Rhw 

e n refr c£< 

Rewriting the latter equation as a function of the potential vector A = £/{2iuj) and using 
for the matrix element the equation 4.2.3 we obtain: 

, . H 2 q 2 f 1 dS\e • M vr \ 2 
aw = 2tt / - — — ' " c| 5.2.36 

e n refr m^ujc J 87H \V k (E c - E v )\ Ec _ Ev=Tiw 



If we assume that the matrix element is constant as a function of k, we can pull it in front of 

1 Aivk 2 rn 
8tt 3 H 2 k 



the integral, and the integral itself gives the value A- ^%T r \ writing the result as a function 



of the reduced electron-hole effective mass m r we have: 

Q 2 /27TI \ 

a(w) = - —i^f) \he-M vc \ 2 ^hw-E G (5.2.37) 

lTieQn re f r mQUJC V fl z J 

In this approximation, the term a(uj)u is proportional to y/hw — Eq, which is the reason it 
is displayed under this form in Fig. 4.7. The Matrix element |e • M vc \ is given by 

h\e-M vc \ = e- f dh^ F ul u {r)pe- i%p u c% {r) (5.2.38) 

J crystal 

The integral can be carried over using the rule for the application of the p operator on the 
Bloch wavefunctions: p — > (hk + p) that yields: 

h\e ■ M vc \ = e • - J— f d 3 re^ul^e-^(hk + p)^(r) 

"crystal J crystal 

= e- y^—{ I dMK^hkuJ + I u^pujf)) 

v crystal J crystal J crystal yO.Z.Oi)) 

Pcv 

assuming e/ /p. The first part of the integral is zero because of the orthogonality of the basis 
Bloch function. Because of the k • p formalism, we know the value of this integral that is 
rather constant for III-V semiconductors as the Kane energy E P = 2p 2 . v /m is about 20eU. 
We may also express it in terms of a dipole matrix element r cv using the relationship between 
the r and p matrix elements (r cv = m ~J Pcv) an d yield 
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To find the accurate values of the matrix element p cv for a given polarization direction, one 
should take advantage from the the fact that the interband matrix element has the same form 
as the one used for the k • p derivation, by just identifying the light polarization direction 
e with the k and using directly the matrix shown in Fig. 5.4. The spin selectivity of the 
|3/2,3/2) to the conduction band 1 1/2, 1/2) is already apparent in the the matrix. 
Using k ■ p approach, we can also explain why the matrix element is a slowly varying function 
of k: as we move away from the k=0 condition, the states will be admixtures of the Bloch 
edge wavefunctions. The respective contribution will scale like the SE/A for the spin-split 
band and SE/Eq for the condution band part. 

Finally, one can rewrite Eq 5.2.37 using the relationship between p cv and the Kane energy 
Ep. To be accurate, on also should consider that with the spin-orbit coupling in effect, | of 
the oscillator strength will be at the first band edge, as shown by the ratio of the contribution 
to the effective masses of the bands. The absorption then writes: 



a(u) = - q - (-rf) -^^hw-E G 5.2.41 

D7re n re f r m c V nr J nuo 



5.3 Effect of carriers 

A key feature of a semiconductor is the possibility it offers to modify its free carrier density 
either by doping, electrical injection or optical pumping. Not only will these carriers modify 
the conductivity of the structure, they also greatly affect the absorption edge. 
When computing the net rate of transition between in a general two level system, the net 
rate of transition writes: 

W net = W abs - Wstim = w(rii - U 2 ) (5.3.42) 

where w is the rate per carrier. We will modify the relations in a similar way for the 
semiconductors. We will assume that the probability of finding an electron at energy E in a 
band is given by a Fermi-Dirac distribution function such that 

^ C exp(^ £ ) + 1 ^ ^ 

In these equations, we choose to allow a different chemical potential for the electrons in the 
valence band and in the conduction band. As the thermalization of the electrons occur on a 
much faster time scale within the bands (typ lOps) than between the bands (typ time Ins), 
this approximation enables us to treat carrier injection in semiconductors. Note that the 
Fermi distribution for the valence band /„ is the one of the electrons, the one for holes being 
simply fi h) = 1 - fi e) . 

As in the two-level system, we write the net absorption as the difference between absorption 
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(proportional to /„(! — f c )) and stimulated emission (proportional to / c (l — /„)), yielding 



2„2 



a{uj) = 2n- 



K 2 q 



dS\e ■ M v 



e n refr m 2 uc J 8n 3 \ V k (E c - E v )\ Ec _ Ev=hjJ 



(/«(l-/c)-/c(l-/«)) 



2tt- 



h 2 q 2 



dS\e ■ M v 



e n re f r mlujc J 8n 3 \W k (E c - E v )\ Ec -e v =Hu> 
which can also be written in the parabolic case: 



(/„(£ M (M) - /c(s fe ,c(M) 



(5.3.45) 



a{uj) 



/2m \ 3 / 2 Ed 

(-^j -£y/tuv-E G (f v (E ktV (tiw)) - f c (E k>c (tuv)) (5.3.46) 



6ii6on re f r moUJC \ h 2 J hw 

where the energy dependence of the Fermi distribution has been explicited in the last equa- 
tion in the kinetic energy terms E ky ( c ,v)- The latter ones can be easily extracted from the 
fact that we are dealing with vertical transitions in k-space: 

E ktV (tw) + E ktC (hu) = hu - E G (5.3.47) 

which for parabolic bands can be expressed as: 

mt 



E k)C 
E k)V 



m* + m* h 



1 (tuv - E G ) 
(hw - E G ) 



(5.3.48) 
(5.3.49) 



5.3.1 Burnstein shift 



In the case of doping, we usually have that if the doping of type n then f v is unity, in the 
case of p doping then f c is zero. The net effect of doping is a shift of the bandgap as the 
band is filled by carrier. 
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Figure 5.8: Burstein shift of InSb (very light electron mass of 0.013 mO) with n doping) 
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This phenomenon has been recently used to fabricate low loss waveguides for InAs-based 
QCLs at A = 3/im, where InAs would normally be strongly absorbing. 
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Figure 5.9: Waveguide including heavily doped InAs layers 

It can be shown, using the Kramers-Konig relations and assuming parabolic bands, that the 
change in absorption edge with the presence of free carriers induces a change in refractive 
index that can be approximated by: 

M*0 = n e , h (5.3.50) 

n re f r m Eq\Eq - E 2 ) 

where Eq ~ Ec + kT is the energy at which the integrand in the KK relation has a maximum. 
The latter relation is valid far enough from the band edge. 

5.3.2 Bandgap shrinkage 

Coulomb interaction introduce a bandgap shrinkage proportional to the n 1 / 3 and is indepen- 
dent of effective mass. This term is difficult to compute accurately but should be taken into 
account for reliable computations of the effective bandgap. 



5.4 Gain 

In the case of non-equilibrium bands (/i c ^ /i v ) then the possibility arises to observe a 
negative absorption, i.e. gain. The condition for the observation is that 



fv - fc < 



(5.4.51) 



5.4. GAIN 



49 



which, when explicitely writing f c and /„, yields that 



which implies that 



exp( ) > 1 



{n c - n v ) > fko. 



(5.4.52) 



(5.4.53) 



The separation between the quasi-Fermi levels must be larger than the photons. (It also 
means that the applied voltage on the diode must be larger than the photons, also!). This 
condition is referred as the Bernard-Durrafourg condition. 

The gain as a function of injected carrier density may be directly computed by using Eq 5.3.46 
and the neutrality condition: 

n = p (5.4.54) 

for undoped material, or more generaly 

N+ + NX+n + p = (5.4.55) 

Using the Fermi distributions and the neutrality condition for undoped material, the dif- 
ference between the separation of the quasi-Fermi levels and the bandgap of GaAs can be 
obtained and is plotted in Fig 5.10. As shown in the latter, for these parameters the trans- 
parency is reached at a carrier density of about 1.1 x 10 18 cm~ 3 . 



Fermi-level difference vs n T=300K 
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Figure 5.10: Separation of the quasi-Fermi levels and the bandgap of GaAs as a function of 
carrier density 

Using the parameters of GaAs {Eq = 1.5 eV, E P = 22eV, and the usual masses) the gain at 
T = 300K has been computed and is shown in Fig. 5.11. 
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gain vs An , T=300K 
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Figure 5.11: Computation of the gain for various carrier densities between 5 x 10 17 to 2 x 10 18 , 
as indicated 



In Fig. 5.12, the energy range where the gain was displayed was broadened to show the energy 
range at which the change of absorption vanishes. The change in absorption with injected 
current, over this wide frequency range, is responsible for the large change in refractive index 
at the gain maximum. 
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Figure 5.12: Computation of the gain for various injected densities, as indicated. The range 
is choosen to show the energy where the absorption is independent of carrier density 
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5.5 Spontaneous emission and luminescence 

To compute the net absorption, we had to deal with the balance between the rate of upward 
transition and downward transitions. Let us write again the net transition rate as: 

Wnet = PBBB cv (f v - f c ) (5.5.56) 

To derive an Einstein-like argument about the relationship between absorption and spon- 
taneous emission, we will be interested in the case that the radiation illuminating p BB the 
system is produced by a blackbody. We will write the transition rate using Fermi's golden 
rule: 

B cv = ^\m nt \j)\ 2 = ^Q 2 £ 2 r 2 cv (5.5.57) 

where we will associate the electric field produced by the single photon in the volume of our 
system, i.e. setting 

1 

We also use the relationship between the interband dipole matrix element r cv and the Kane 
energy E P 

ri = 7^1? , (5-5.59) 



^ = o n re/reo£ 2 V. (5.5.58) 



B cv = o -w-. (5.5.60) 



as a result, we obtain 

itq 2 h Ei 
2n 2 efr € m E G ' 

The blackbody electromagnetic energy density, usually expressed as a function of frequency, 
must be expressed in density of photon per unit energy, i.e. 

8nn 3 refr h 2 u 2 V 1 

PBB ~ exp(-&)-l- (5 - 5 - 61) 

We then write a detailed balance between stimulated emission, absorption and spontaneous 
emission at thermal equilbrium: 

PbbB cv J v (\ - f c ) = pBBB cv f c (l - f v ) + A cv f c (l - f v ) (5.5.62) 

Using the definitions of f c , f v as above, we obtain that A cv and B cv must be related by 

A cv = B cv ^ . (5.5.63) 

The coefficient A cv has to be interpreted as the rate at which a carrier in the band, having a 
hole to recombine with, will incur spontaneous emission. We define the inverse of A~^ = t$ 
as the spontaneous emission lifetime Ts- Finally, using the value of B cv derived from above, we 
have: 

1 = 9 f 2 "T /r E G E P . (5.5.64) 

A typical value for the spontaneous lifetime is t s = 0.7ns for GaAs. As shown by the above 
equation, this time will drop as the energy of the gap increases. 
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5.5.1 Relationship between absorption and luminescence 

Note that we can write the absorption coefficient a as: 

a = ^B cv (f v -f c ). (5.5.65) 



The Einstein relation between A cv and B cv allow us to write a very general relationship 
between the spontaneous emission and the absorption. The spontaneous emission rate at 
the photon energy hv 

A cv fc(]- — fv) 



spon 



B C v ( fv fc 



a(hv) 



aihv) 



c 8nn 3 ref y / c (l - f v ) 

fv fc 

1 



n re f r 



he 2 



(5.5.66) 



lim 2 refr h 2 v 2 



hh 



'c* exp(-j^) - 1 

For large gap materials, and at reasonably low temperatures, E G » kT and therefore 



spon 



a(hv) 



^n 2 refr h 2 u 2 



h 3 c 2 



exp{- 



kT' 



(5.5.67) 



5.5.2 applications: luminescence lineshape 



In general, we saw that for Bulk materials, the absorption has the form \JTiuj — Eq with 
some weaker energy dependences. As a consequence of Eq. 5.5.67 the luminescence lineshape 
should conserve the \[Wo — Eq behavior at the low energy side and an exponential decrease 
at the high energy side. 

C(hv) w ^hv-E G ex V (-^) (5.5.68) 



The width of the curve should be approximatly 1.8kT. 



Measurement 



Photoluminescence of GaA5 at 300K 
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Figure 5.13: High energy tail of the photoluminescence of GaAs as a function of pumping 
intensity. The fitted temperature is indicated. 
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The photoluminescence of a undoped GaAs at room temperature is shown in Fig. 5.13 as 
a function of wavelength, along with a fit that uses an expression equal to Eq. 5.5.68 but 
expressed in wavelength units. 



5.5.3 Electronic temperature 

The measurement of the high energy tail of the distribution, because of the relative slow 
variation of the absorption, enables the measurement of the carrier temperature. This was 
recognized very early by Jagdeep Shah of Bell laboratories who first measured the shape of 
the luminescence as a function of optical pumping intensity. Result of such an experiment 
is shown in Fig. 5.14 




Figure 5.14: High energy tail of the photoluminescence of GaAs as a function of pumping 
intensity. The fitted temperature is indicated. 



A plot of the carrier temperature as a function of pumping intensity demonstrated the key 
role of the optical phonon energy in the carrier energy dissipation: 



CHAPTER 5. BULK SEMICONDUCTORS: BANDSTRUCTURE AND FUNDAMENTAL 
54 GAP 




Figure 5.15: Electron temperature as a function of pumping flux. The activation energy 
found (33meV) is very close to the optical phonon energy (36meV) 



This technique became a standard procedure for the measurement of the carrier temperature 
in all optoelectronic devices (lasers, transistors, etc..) 



5.5.4 Gain measurement 

An other interesting application of the relationship between luminescence and absorption is 
the possibility of measuring gain from the measurement of unamplified luminesence. In fact, 
assuming that f c and /„ are characterized by quasi- Fermi chemical potentials /i c and fi v , the 
ratio / c (l — f v )/(fv ~ fc) appearing in Equ. 5.5.66 does yield: 



(5.5.69) 



fc(l ~ U) _ 1 
fv-fc " exp( ^fa^ ) - 1 



in which case the absorption is related to the luminescence by: 



«(H = r »°» 8im * refr hw exp( kf } " 1 (5 - 5 - 70) 



and the latter equation can be used to measure the absorption (and therefore the gain) 
through a measurement of the unamplified luminescence as a function of injected current. 
This technique was first demonstrated by C. Henry from Bell Labs. 
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Figure 5.16: Schematic drawing of the luminescence measurement geometry- 




Figure 5.17: Luminescence as a function of injected current. The measurement of the gain 
is deduced from the ratio of these curves 
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Figure 5.18: Measured gain at 1mA and 64mA injection current 
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Figure 5.19: Gain at the laser wavelength versus carrier density 



5.5.5 Bimolecular recombination 

Let us consider the term / c (l — /„) at low injection levels, in the classical regime, we have 



fc 



1 - fv = 1 



exp(^) + 1 



cxp( 



exp(^=^) + 1 
As a result, the product / c (l — /„) is given by 



E c — A* c x 
kT } 



(5.5.71) 



(5.5.72) 
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We then want to compute the total recombination rate from all the carrier population. To 
this end, we integrate the spontaneous emission as a function of photon energy: 

R S pon = I —pj(hu)f c (l - f v )dhv 

J T~S 

exp(^) 1 /2m r x3/ 2 [<*> r- hv ^ 73) 



Ts 2tt 2 



exp(^y^) 1 /2m r A;T\3/2 



t s 2tt 2 



y — ^ — J J VMexp(-M)dM. 



The integral J °° y / wexp(— u)du = The above equation may be simplified by expressing 
the quasi-Fermi level implicitly through the electron and hole concentrations, always in the 
classical limit: 

E c - fi c 



n = iV c exp( 

,E V - fi v , 
p = N v exp( kT ) 

where the quantum concentration N c>v is given by: 



(5.5.74) 



The product np yields (law of mass action) 



np = N c N v e^(- {Ec - K) - T ^-^ ) (5.5.76) 
That allows finally us to write the total spontaneous emission rate: 

Rspou = 1 ^ P N N j (5.5.77) 

written usually as a function of the bimolecular recombination coefficient B 

R S pon = Bnp (5.5.78) 

As the value of B is given by: 

B = WW7- (5 ' 579) 

c v spon 

B is both material and temperature dependent. The temperature dependence is T~ 3 / 2 and 
the material dependence can be expressed using the known values of the quantum concen- 
trations and of the spontaeous lifetime. 

It should be remembered that the range of concentration in which this expression is valid: 
it assumes classical distribution (limitation towards the high concentrations) but does not 
include the excitonic effects and so is expected to fail at very low concentrations. 
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Chapter 6 

Quantum wells and nanostructures 



6.1 Envelope function approximation 

The problem we will try and solve now is the one of a heterosructure, in which two materials 
A and B are sandwiched together. Of course such a material could be seen as a new material 
by itself, and its band solved by ab initio techniques, but such a computation is very heavy, 
time consuming and moreover does not give much physical insight into the result. The 
envelope function approximation solves this problem in a very efficient and elegant manner. 
It is widely used to predict the optical, electrical properties of semiconductor nanostructures. 

6.1.1 Multiband case 

At the core of the envelope function approximation is a generalization of the k • p approx- 
imation: it is postulated that the wavefunction can be written as a sum of slowly varying 
envelop functions f^' B (r) that will modulate the Bloch function of the material, namely: 

*(0 = £/i* B <£( r )- (6-1-1) 
i 

Behind the equation 6.1.1 is the idea that at each point, the wavefunction is described by a 
k • p decomposition and that this decomposition depends on the position. Furthermore, it 
is assumed that 

1. the envelop function ff ,B (r) is slowly varying compared to the Bloch wavefunction, if 
ff' B {r) is written in a Fourier decomposition, the wavevectors are close to the center 
of the Brillouin zone 

2. the Bloch functions are identical in both materials, i.e. u^ k (r) = u^ ko (r). This also 
implies that the interband matrix element (S\p x \X) is equal in both materials. 

It allows us to write te wavefunction as 

nr) = J2tf' Bu ^(r)- (6.1.2) 
i 
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Let us assume first a quantum well, in which a layer of material A is sandwiched into a 
barrier material B. Because of the in-plane translational invariance, the wavfuction may be 
written as plane waves: 

fi(r±,z) = -j=e*^ X i{z) (6-1.3) 

where z is choosen as the growth direction and k± = (k x ,k y ) is the in-plane wavevector. 
Note that this convention there is confusion in the litterature as the sign _L may either mean 
perpendicular to the plane of the layers or to the growth axis. The Hamiltonian is then 

H = ^ + V A (r)Y A + V B (r)Y B (6.1.4) 

where the functions Ya(z) and Y B (z) "turn on" the potential in the respective layers. We 
will develop our system close to k=0. To solve the system, we must: 

1. Let H act upon \l/(r) 

2. multiply on the left by u^ n0 (r)e~ tk±r± Xm( z ) 

3. integrate over space 

We have to use the following relations. As the envelop function is slowly varying, we may 
write 

/ d 3 rftf m u^u m = fif m / d 3 ru^u m = ft f m 5 lm (6.1.5) 

J cell J cell 

and take advantage of the fact that the band edge are eigen function of the Hamiltonian at 
(k=0): 

+ V A ' B (r))u mfi (r) = e£X,o(r). (6.1.6) 

The derivation is rather tedious, but one should note the similarity with the normal k ■ p 
technique by considering the action of the operator p on the wavefunction: 

p(e k ^ X i(z>i(r)) = (hk ± - ih-^+p)e ik ^ X i(z)Mr) (6.1.7) 
and we then may consider the substitution 

p -> (hk ± - ih— + p) (6.1.8) 

oz 

where it is understood that p then acts only on the Bloch part of the wavefunction. Using 
the above substitution into the Hamiltonian, one finally get the following set of diferential 
equation written in a matrix form: 

D(z,-th-^)x = e X (6.1.9) 
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where the elements of the matrix D_ are given by the equation 

, fc2p fc2 Q2 



5, 



m 



2m 2m <9z 2 

fi/c_L ... . . ih ..... d 
f _(Z|p ± | ro >-_(Z|p»_ 



(6.1.10) 



and the matrix element (Z|p|m) = § u u*pu m d?r. 
6.1.2 One Band model 

As an example, let us consider first a pure one band model. Setting 

efY A + efY B = V(z) (6.1.11) 
we obtain the Schrodinger equation: 

( -is +/w w ,1= « (, i (6 ' L12) 

This equation is of course the one of a free electron in a potential. As a consequence, it does 
not take into account the dispersion of the band. A much better model would be to replace 
mo in equation 6.1.12 by the band effective mass m e //. For an isolated band such as the 
heavy hole band, it is a rather good approximation. Formally, it should be done by adding 
the effect of the remote bands in the matrix 6.1.10, the result of which is given (see Bastard) 
for a 8 band model with the remote bands as 



/ h 2 k 2 h 2 d 2 Hk\ ih d 

V ( [4 + V m (z) + — - 75— Tj-aRm + (1\P±H (l\Pz\m) — 

^ V L 2m 2m az 2 J m m cte 

m=l 



2dzM?*dz 2 ^ L Mf%dz dzM™ J 2 ^ Mr 



where effective mass parameters are defined as 



(6.1.13) 



One sees that for a one band model, neglecting the in-plane dispersion, the term %-§-irkz-§- 
is the only one remaining and will then change the effective mass. Boundary conditions can 
be derived and be shown to force the continuity of the wavefunction and of the quantity 

—fysT (6-1-15) 

m(z) oz 



which is proportional to the probability current. 
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Using this model, on immediatly obtain the confinement energy in quantum wells, that will 
express themselves by a increase of the apparent conduction band and valence band energy. 
This effect will be stronger for electrons that are lighter than for the heavier holes. 



\ 1 \ 1 \ 
IH \P= 2 \P= 3 




50 100 150 

l (h 



Figure 6.1: Computed confinement energy of electrons in the conduction band as a function 
of well width. 




Figure 6.2: Computed confinement energy of the valence band as a function of well width. 



6.1.3 Two band model 



For the conduction band, a very nice model is the one in which one keeps one valence band, 
creating a two-band model. For simplicity, let us look at the states at kj_ = and neglect 
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the free electron term. We obtain the system of equation given by: 

ih d 

Pea ~7T~ 

mo dz 



V c (z)f c PcvT^fv = e/c (6.1.16) 



1 P* cvir fc + V v (z)f v = ef v (6.1.17) 



m cv dz 

Extracting f v from the second equation yields: 

= vm=i^tJ< (6118) 

replacing into the first equation, after substitution, the following reslult is obtained. 

] ' 1[,kX ' () ' 9 f c + V c (z)f c = ef c . (6.1.19) 



m dz e — V v (z) dz 

2 2 

uciiiiiLiuii ui me ivaiic cncigj -t^p — ; " 

effective mass 



Recalling the definition of the Kane energy Ep = ^Pc V an d defining an energy-dependent 



m(e,z) m e — V v (z) 
we obtain finally a Schrodinger-like equation 

h 2 d 1 d 



1 1 E P . 

(6.1.20) 



2 dz m(e, z) dz 



fc + V c (z)f c = ef c . (6.1.21) 



This model is very useful to model the electronic states in the conduction band with the 
inclusion of the non-parabolicity. It is very widely used in the study of intersubband transi- 
tions. 




Figure 6.3: Energy states of a quantum well computed with a two-band model, and compared 
with a one-band model (dashed lines). The growing importance of non-parabolicity as one 
moves away from the gap is clearly apparent. 



64 



CHAPTER 6. QUANTUM WELLS AND NANOSTRUCTURES 



6.1.4 Full model: the valence band 

To be accurate, the full kp envelope function model must be solved, that is the set of 8 
differential equations defined by the equation 6.1.9. Solving the full model is essential in the 
valence band because no simplified model can easily be used if some degree of accuracy is 
sought. It is observed that: 

• It is usualy convienient to use the growth direction as the quantization direction for 
the angular momentum. 

• The confinement potential lifts the degeneracy between the heavy and light hole bands, 
because their different mass induce a different confinement energy 

• The in-plane dispersion is highly non-parabolic because of the coupling between the 
bands induced by the in-plane momentum. In particular, this coupling prevents any 
crossing between the light-hole derived and the heavy hole-derived band (see dahsed 
lines) 

• In some cases, the mass is inverted: the bottom of the LH1 band has a electron- 
like character over some portion of reciprocal space because of the repulsion and its 
proximity to the HH2 state. 

This effect is shown schematically in Fig. 6.4. As a result, the computed band structure is 
usually fairly complex, and yield results such as the ones shown in Fig. 6.5 



Bulk (k=0) 



J - % 



j,« t'/i. ■ h 
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En ■ 
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J 2 ~- t % 
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No interaction 




^(fi 1 /2mJ(-r,-i ! ) 



(c) 



With interaction 




(d) 



Figure 6.4: Schematic description of the origin of the valence band dispersion in the quantum 
well showing schematically the effects of confinement and interactions. 
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Figure 6.5: Dispersion of the state of a quantum well in the valence band. 

6.2 Absorption: interband case 

The states computed using the envelope function approximation can now be used as a basis 
for the computation of the absorption in the system. Let us use here the dipole interaction 
Hamiltonian (the reason for using it will appear later) so that 

H int = -qS-r (6.2.22) 

and let us consider the transition from a valence band state given by 

\nk ± ) = -J- Xn ( z )e ik ^u v (r) (6.2.23) 
VA 



to a conduction band 

\mk' ± ) = -^= X m(z)e ik '^u c (r) (6.2.24) 



where r = (r±, z) is the position and r_|_ is the in-plane coordinate and A the sample's area. 
The matrix element is then 

(mk' ± \£ ■ r\nk ± ) = ^ I d 3 r X * m (z)e- ik '± r u* c {r)£ ■ r X n{z)e iqr e ik ^ r u v {r). (6.2.25) 

J space 

The above integral can be simplified by: 

• Breaking the integral into a sum of elementary cells ^ i i over which the integral is 
carried over 

• Assuming that the envelope functions are slowly varying and therefore can be pulled 
out in front of the integrals. 
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We then obtain: 

(mk' ± \£ ■ r\nk ± ) = 1 £ xU^Xn^)^- £ e^~ k '^ 



-M/ d 3 r'u:(r')S-r'u v (r)e^ +<1 - k ^ r ' (6-2.26) 

Vcell \- J cell 

+£-r' I d 3 r'u* c (r')u v (r)e i ^ +q - k ^ r ' . 

J cell 

that can then further be simplified (note that the second integral is zero because the Bloch 
functions are orthogonal): 



/oo 

X*m,cXn,v (6.2.27) 
■oo 

The corresponding scattering rate, computed using Fermi's golden rule, writes: 

W :q 2 \r co \ 2 \(m,c\n,v)\ 2 . (6.2.28) 



As in the three dimentional case, it should be summed in k-space for all available states that 
satisfy the k-selection rule, and multiplied by the corresponding Fermi distributions. We 
obtain for one pair of subband: 

w v ^ c = —q 2 \r cv \ 2 \(m,c\n,v)\ 2 -^Q(hw - (E G + e c , m + e„ >n ))/ W)n (l - / c , n ) (6.2.29) 

Notes on this result: 

• The scalar product | (m, c\n, v)\ 2 will be zero for quantum well states of opposite parity 
in symmetric quantum wells. 

• In general, this term will be very small for states with different level index. 

• The matrix element r cv depends on the polarisation direction and on the nature of the 
valence band Bloch states. 

Let us consider the last point a little more in detail. For example, let us assume we are 
looking at a transition between a state ||,|) and the conduction band \i§ |). Using the 
relation between the dipole matrix element and the momentum matrix element: 

r cv = — % ^rVcv (6.2.30) 
m E G 

we can then use the Kane matrix k ■ p 5.4 to obtain the value of the matrix element, that 
would be -^Pcv for this heavy-hole to conduction band transition in the TE polarization 
direction. One can check that considering a light hole to conduction band transition would 
lead to a coefficient equatl to 
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Table II. — Selection rules for interband (inter-subband) transitions obtained from the 
absolute value of the matrix element (S t \z.p \ u a ) , a = T g , r 7 .k L = 0. 



Polarization 








Type of transitions 


Propagation 
parallel to z 


n 
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Figure 6.6: Interband selection rules for quantum wells 
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Figure 6.7: Quantum well absorption in TE and in TM, showing the absence of HH transition 
in TM 

Finally, the absorption can be obtained from the scattering rate by computing the ratio of 
the absorbed power divided by the incident one. For a quantum well at an incident angle of 
9, the result yield: 

hww v ^ c q 2 m r huj 2 2 , 

a = ~ n 7m = is 7^\Vq- r cv\ \{m,c\n,v}\ (6.2.31) 

±e nc£ z cos(V) n 6 e n re f r c cos(fl) 

assuming the whole beam is intercepted by the quantum well. In the above equation, the 
absorption is number, and represents the fraction of the beam absorbed by the quantum 
well system. For situation where the incidence angle is 9 = it/ 2, the latter assumption is not 
possible and one then defines an absorption coefficient given by: 

=t3 \Vq- r cv\ \{m,c\n,v)\ — (6.2.32) 
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where d qw is the thickness of the quantum well (but can be defined in an other manner as 
the period length for a superlattice) and the overlap factor of the quantum well with the 
intensity of a guided modeF" is defined by: 



f-t% E2dz 



J — oo 



(6.2.33) 



Note that the product d qw T is almost indepdent of d qw as long as the quantum well is much 
smaller than the width of the guided mode. 

6.2.1 Interband spectroscopy: notes on the techniques 

Comparison between the different techniques used. 

• Absorption 

• Photoluminescence 

• Photoluminesence excitation 




1.70 1,90 
Photon energy [eV| 



Fig. 7.16a-t. Comparison between 
absorption and PLE spectra in 
GaAs/AlGaAs quantum well (QW; 
samples measured at low temper- 
ature, (a) Absorption spectra in 
a 20 nm well showing structures 
due to excitonic absorption peaks 
(labeled n = 1.2, etc.) [7.45]. 
(b) PLE spectra from a sample 
containing two QWs with widths 
of 5 and 10 nm (shown in lower 
and upper panels respectively). 
The broken curves show the PL 
from these two wells [7.46]. (c) 
PLE spectra from a p-doped well. 
The broken curve was obtained 
by selling the spectrometer to ad- 
mit only excitonic emission (indi- 
cated by the broken vertical ar- 
row). The vertical arrow labeled 
Eg marks the energy gap of the 
AlGaAs harrier layer. The solid 
curve is the PLE spectrum ob- 
tained when the spectrometer was 
set at the energy of the electron- 
to-neutral-acceptor emission (in- 
dicated by the solid vertical ar- 
row) [7.47 1 



Figure 6.8: Comparison between absorption, PL and PL excitation (From Cardona's book) 
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Figure 6.9: Comparison between absorption, PL and PL excitation (From Weisbusch and 
Winter) 



6.2.2 Note on the QW absorption 

Let us assume the transition between a pair of subband in a quantum well. The absorption 
for normal incidence 9 = it/ 2 (considering only the heavy hole transitions first): 

[rS ,|2| ( i )[a (6 . 2 . 34) 

We should first note that the overlap factor | (m, c\n, v ) | 2 ~ 1 for a transition between two 
ground states. Then, we shall use the relationship between the dipole matrix element and 
the momentum matrix element 

— ih 



m E f 



-Pcv (6.2.35) 



G 
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and we have, for a transition between hhl and El: 

For the same transition, we may also express the joint density of state effective mass as: 

l r,hh 



m r,hh = m; 1 + mil (6-2-37) 



For the heavy hole state, we will use the in-plane mass of the ||, |) state that gives (see the 
kp matrix Fig. 5.4) as 



m hh = ~ m o 1 + ~et = ~ m o 1 + m o 1 7rf- (6.2.38) 



For the electron effective mass, we have: 



-i — i '-/-.'/■ f E P 



m; 1 = m Q 1 + f_ + —f _ (6.2.39) 

3m E G 3m (E G + A) v ; 



that then yields for the reduced mass: 



m r,hh = m - 



i Ep \ 3Ep 

™o Up! + wTp! _l a J ~ o^TWZ (6.2.40) 



-6 E G Smo^c + A) 7 2m 0J E G 



Substituting both reduced mass and matrix element by the above expressions in the absorp- 
tion yields: 

otfih = % = " (6-2.41) 

3 x Ahe n re f r c 3n re f r iraf ine 

where a.f ine ~ 1/137 is the fine structure constant. For a refractive index of n re f r = 3.6, 
the absorption corresponds to a hh = 4.26 x 10~ 3 . The effect of the refractive index is only 
there when measuring thick samples; for a very thin quantum well suspended in vaccum the 
absorption would simply be equal to 27ro;/j ne /3. It is interesting to note that not only this 
results does not depend on the detail of the quantum well such as its thickness or barrier 
height; it does not even depend on the value of the Kane energy Ep. 

In most samples, the transition between LH1 and El is very close to the HH1 El transition 
and, because of excitonic effects, it is not possible to find a plateau of absorption between 
these two values. For this transition, the light hole mass (assuming the ||, \) basis state) : 



m ih 



1 = -mo 1 + ^- = -mo 1 + mo 1 ^- (6.2.42) 

The reduced mass is, in turn: 

-l -i/5-Ep Ep s 5E P 

m,, = m n ; —— ) ~ (6.2.43) 
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Similarely, the matrix element r cv is given by 



- i£fe (6 - 2 - 44) 



Again, combining the two, we obtain: 



a lh = —i = \ (6.2.45) 

I0ne n re f r c o n re f r naf ine 



The sum of the absorption caused by the lh and the hh state is then 



16 1 

a hh ,i h = — (6.2.46) 

10 Tl re j' r 7lOi fi ne 



which is equal to 6.8 x 1CT 3 for a refractive index of n re f r = 3.6. The data in the litterature 
yield value somewhat larger. As an example, the data from Stolz et al for various well width 
in InGaAs/AlInAs samples yield a value of (0.82-0.89%), and an analysis of the absorption 
in InGaAs/InP quantum wells of Sugawara yields 0.85%. A enhancement given by the 
Sommerfeld factor, equal to 2 close to the bandgap, is perhaps the reason for the discrepancy. 
If it is true, the value for very low mass materials such as InAs should yield lower values. 
Indeed, values for the GaAs-based quantum wells reported by Masselink, with heavier masses 
and where excitonic effects ar more prevalents, are closer to 1% absorption. 



Additionally, measurement of such small absorption factor are always difficult experiments. 
In particular, baseline issues become critical. The presence of standing wave brought about 
by the reflection on the surface will equally skew the results. 
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Figure 6.10: Sommerfeld factor for the exciton: 3D versus 2D case (Weisbuch and Winter). 
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Figure 6.11: Absorption for quantum wells of various thicknesses 



6.2.3 Effect of strain 



The use of strained layers has been a very important milestone in the development of quan- 
tum well lasers. 
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Figure 6.12: Effect of the material strain on the band structure. The strain anisotropy will 
split the valence band degeneracy between the HH and LH states by a value S 



In the early nineties, it was realized that even thought lattice matching had to be preserved 
at all costs for the growth of quantum well lasers when thick layers were involved, one could 
grow very thin layers pseudo-morphically, i.e. in such a way that the grown material adapts 
its lattice to the one of the substrate. The result is a material in which very large amount of 
strain can be incorporated, corresponding to a lattice mismatch as high as 2%. The resulting 
strain is very anisotropic, as the material is constrained in the plane of the layer but can 
adapt its lattice spacing along the growth axis. The effect of this anisotropic strain is to 
split the valence band maximum and split the heavy hole and light hole bands by a value 
S, as indicated on figure 6.12. As the HH states derived from the ||, |) are heavy along 
the growth axis but light in the plane, the effective density of state in the plane is strongly 
reduced. 
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Figure 6.13: Valence band dispersion and density of state for a unstrained 8nm thick GaAs 
quantum welll 




Figure 6.14: Valence band dispersion and density of state for a strained 8nm thick InGaAs 
quantum welll with 20% Indium 



The figures 6.13 and 6.14, showing the dispersion along with the effective density of state, 
show the significant reduction in the density of state of the valence band that can be obtained 
by about 1% strain in the quantum well. As a result, the density of carrier needed to reach 
transparency, and as a result the transparency current density, are very strongly reduced, as 
shown in Fig. 6.15, 6.16, 6.17. 
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Figure 6.15: Computation of the gain for various carrier densities (indicated in units of 
10 i8 cm -sv bulk Q a As. 
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Figure 6.16: Computation of the gain for various carrier densities (indicated in units of 
10 18 cm~ 3 ): unstrained GaAs QW. 
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Figure 6.17: Computation of the gain for various carrier densities (indicated in units of 
10 18 cm -3 ): a strained 8nm thick InGaAs quantum welll with 20% Indium 



6.3 Intersubband absorption 

6.3.1 First approach: dipole in a one-band model 

We use the same approach as above, using again the dipole interaction Hamiltonian. We 
now consider a transition between two conduction band states: 



to a conduction band 



\nk±) 



1 



Xn(z)e ik±r ^u c r 



Xm( z y k ' ±r± u c r 



(6.3.47) 



(6.3.48) 



where, again r = (r±, z) is the position and r_|_ is the in-plane coordinate and A the sample's 
area. The matrix element is then 



{mk' ± \£ ■ r\nk±) 



A 



d 3 r X * m ( z )e~ lk ' ±r u* c (r)S ■ r Xn (z)e iqr e lk ^u c (r). (6.3.49) 



iqr ik ± r 



space 



As above, we procceed the same way as for the interband case and we: 



Break the integral into a sum of elementary cells £\ i i over which the integral is 
carried over 



Assume that the envelope functions are slowly varying and therefore can be pulled out 
in front of the integrals. 
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We then obtain: 

(mk' ± \S ■ r\nk ± ) ^E^lj^E^"^ 

-M/ dh'u* c {r')£ -r'uciry^-W (6-3.50) 

Vcell Ucell 

+£-r' I d 3 r'u* c (r')u c (r)e l(k ^ +q - k '^ r ' 

J cell 

In this case, it is the first integral that vanishes and the second part that remains: 

(mk' ± \£ ■ r\nk±) = 5(k± + q- k' ± ) ^ Xm( z )Xn(zi)a z £ ■ r (6.3.51) 

iz 

that can be converted into an integral that finally yields: 

/oo 
dzx* m {z)zxn{z) (6.3.52) 
-oo 

As in the interband case, the matrix element should be inserted in Fermi's golden rule; 

w v ^ c = ^q 2 \(m,c\z\n,v}\ 2 . (6.3.53) 

As in the interband case, it should be summed in k-space for all available states that satisfy 
the k-selection rule. Since the subbands are paralell in this one-band model, absorption 
exists only when the pair of subband is separated by exactly the photon energy: 

a(hw) = —q 2 \(m, c\z\n, v)\ 2 .5(E m - E n - hw)(n n - n m ) (6.3.54) 

6.3.2 Absorption in a quantum well: a two-band model 

The problem arising with the previous model is that it does not allow the introduction of non- 
parabolicity. It is valid then only for confinement energies much smaller than the bandgap 
(E n) E m « E G ). One could be tempted to use the same approach when computing the 
absorption in the multiband case, taking as the matrix element the envelope function of the 
electron. This approach was used in some of the works, but lead to difficulties since, strickly 
speaking, the envelope functions of the conduction bands are no longer orthogonal to each 
other ((xm\Xn) <^nm)-A much better approach is to compute the matrix element directly 
in the multiband model. If one wants to treat conduction band state, an effective two-band 
model is perfectly valid. Let us then assume a two-component wavefunction: 

^i = / c (1) « c + /X (6.3.55) 
^>2 = f { c 2) u c + f 2 v u v (6.3.56) 



Let us then evaluate the matrix element p z : 

(*lb,|*2> = Ui 1] U c + fi l) U v \p z \f { c 2) U c + fi 2 \ v ) 

= (f^u c \ Pz \fPu c ) + (f^u c \p z \f^u v ) + (fi l) u v \p z \f^u c ) + {fl l) u v \p z \fi 2 \ v ). 



(6.3.57) 
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Let us consider the first term, it yields: 

(f^u c \p z \fPu c ) = (fP\ Pz \fW)(u c \u c ) + (fP\fP)(u c \ Pz \u c ) 

= (fv\p z \fp) 

Similarely, the second term will yield 

(/iVb z |/iV) = (/ c (1) l/ c (2) )p c , 



(6.3.58) 



(6.3.59) 



Using similar derivations for the third and fourth terms, the result can be summarized in a 
matrix form: 

Pz Pcv 
Pcv Pz 



p 



(6.3.60) 



acting on the components (f c , f v ) of the wavefunction. Dropping the diagonal terms as 
Pz « Pcv, we finally obtain: 



(*ib*l* 2 > = (/ c (1) l/i 2) )^ + (fP\fP)^ 



(6.3.61) 



Using the relationship between f c and f v given by the equation 6.1.18 and using the definition 
of the energy- dependent effective mass 6.1.20, we finally obtain: 



(* 1 bJ*2> = (/«iH^) '"" 



dz J m(E,z) m{E,z) { dz )lJc 1 



(6.3.62) 



In this picture, both intersubband and interband transitions are treated on the same footing. 



6.3.3 Experimental results 

As mentionned above, the dipole matrix element for intersubband transition is non-zero 
only for the z-component of the electric field. As a result, a number of geometries have been 
developed to measure the absorption, as shown in Fig. 6.18. 



Multi-pass A 



Single-pass A 



Mesa 



Multi-pass B 



Single-pass B 




Brewster 



Figure 6.18: Experimental geometries allowing the measurements of intersubband transitions 
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An exemple of absorption between two subband is shown in Fig. 6.19. 




PHOTON ENERGY (meV) 



Figure 6.19: Intersubband absorption between two bound states 

In contrast to interband transition, where the interband matrix element r cv is the dominant 
term, the atomic-like nature of the joint density of state as well as the tailorability of the 
potential enables the fabrication of complex energy ladder structures. 
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Figure 6.20: Intersubband absorption in a multiquantum well designed for triply resonant 
non-linear susceptibility 
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When the excited state is in a continuum, the absorption is broadened. 




Such quantum well absorption is the basis for the so-called Quantum Well Infrared Photo- 
conductors (QWIP) devices. 



It is even possible to create quasi-bound states in the continuum using Bragg reflection, as 
shown in Fig. 6.22. 
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Figure 6.22: Bound-to-quasibound transition. The confinement of the upper state is created 
by Bragg electronic Bragg reflection 



6.3.4 Intersubband absorption: multiband problem 



Intersubband transitions are also possible in the valence band. However, the situation is 
much more complex because transitions with an interband character (i.e. where the Bloch 
part of the wavefunction is changed, such as a transition between HH1 and LH1 states) will 
coexist with transitions with a more intersubband character (like a transition between HH1 
and HH2). The matrix element is in general dependent on the in-plane wavevector, further 
complicating the analysis. However, a multiband k ■ p analysis yields spectra that fit very 
well the epxerimental data. 
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Figure 6.23: Intersubband absorption measurements for both light polarization directions 
TE and TM in a SiGe quantum well. Due to the presence of states with different Bloch 
wavef unctions, the absorption is allowed for both polarization directions, a) Experimental 
result, b) Computed absorption using a 6 bands k ■ p approach, c) Band structure of the 
quantum well, indicated the location of the different confined states as well as their main 
character (HH, LH, SO) 



Chapter 7 

Electric and Magnetic Fields 



7.1 Franz- Keldish effect 



The electric field S is applied on a bulk semiconductor of gap Eq. If we want to consider 
the most general situation the complete Schroedinger equation must be solved (not taking 
into account excitonic contribution): 



h 2 

V 2 - eFz ) Vn(r) = E n ^ n {v) (7.1.1) 



2m r 



The solution of this equation can be expressed as a combination of Airy functions. The 
oscillatory character of the Airy functions introduces oscillations in the absorption above 
bandgap. 

The low energy absorption edge is modified according to the following formula: 



otFK — ol{uj, S) ~ 



8tt> c 



£ 



exp 



4 

3F 



'2m r (E c 



h 2 



My 2 ' 



(7.1.2) 



Absorption now extends below E g and the profile is modified (see continuous line Fig. 7.1 ) 
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Figure 7.1: Absortpion edge calculated for bulk (continuous line) and quantum well for an 
applied electric field of 100 kV/cm 

If excitonic effects are introduced, the absorption profile is modified according to the following 
formula: 

a(u) = a FK ^r(l - -^=)exp {^ ln (~^~ ))) (7.1.3) 

The main effect on the excitonic resonance is the blurring of the excitonic absorption peak 
as a function of the applied electric field. 



7.2 Quantum confined Stark effect 

The quantum confined Stark effect arises from the application of a DC electric field on a 
quantum well system. Two main configuration are analyzed: the electric field is parallel to 
the growth axis of the well or it is perpendicular to the growth axis. 

7.2.1 Transverse QCSE 
Interband case 

When the applied electric field is sufficiently low, a perturbative approach can be followed. 
Let us evaluate the effect of an electric field F on a quantum well of width L. We consider 
the ground state in the conduction band and the first heavy hole state in the valence band. 
If we label the states as , |^) we can write the correction to the energy levels due to 
the perturbation W = —eFz. At the first order the energy correction to the levels will be 
zero because the quantum well is symmetric and the centroids of the wavefunctions will have 
the same expectation value: 
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ae 1 ^ = mwm - mwm = esm\zm - m*wi)) = o (7.2.4) 

We then consider the second order correction El 1 — V" 1 . , KM^fei)] ■ nly the nearest levels 

n L^ij^n E n -Ei J 

will be included in the calculation and the energy correction will then be: 



AE ¥tark = -e^ 2 



m\zM)\ 2 , m\z\w)\ 



+ (7.2.5) 
E% - E{ E v 2 - El 

The application of an electric field results in a redshift of the transition: the effective bandgap 
of the quantum well is reduced. If we employ the infinite well approximation where the energy 
levels are expressed by E n = n 2 „ r f n2 r2 and the matrix element is Z21 = tt-t we can express 

Zm e,hh Ll y *" 

this energy shift in terms of: 

U%Lk = {j) 2 \{K + rnl h )L A (7.2.6) 

where the strong dependence from the well width is explicit. 

The excitonic resonance is not blurred but can be shifted for high values of the applied 
electric fields (see Fig. 7.3): this effect is at the basis of the electroabsorption modulators. 

Intersubband case 

In the intersubband case, we will obviously have zero shift at first order because of the 
symmetric profile of the quantum well. At the second order the calculation is analogous to 
what seen in the previous paragraph, by taking into account only the ground state and the 
first excited state of the quantum well. The shift in the intersubband case reads (evaluated 
for the conduction band): 

AE StarkISB - _ - - ™ e (7.2.7) 

Note that in this case the shift is positive: the ground state of the QW is lowered and the first 
excited state is lifted up. The intersubband transition in a symmetric quantum well will blue 
shift as a function of the applied electric field F. The range of validity of the perturbative 
treatment is given by AEqcse *C — E° ~ 2f ^ 2 ■ For a typical value of 100 kV/cm for the 
applied electric field, we obtain A ^° c if = 0.1 for a L=12 nm wide quantum well. First order 
QCSE can be observed in asymmetric quantum wells or in coupled quantum well systems, 
like the quantum cascade laser 

7.2.2 Longitudinal QCSE 

In the longitudinal QCSE the main effect is the blurring of the excitonic resonance as ob- 
served in the bulk: this is expected since in the plane of the quantum well the carriers are 
free to move as in the bulk. 
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7.2.3 Quantum Well exciton resonance in an electric field 
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Figure 7.2: GaAs/AlGaAs multiquantum well structure 
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Figure 7.3: Absorption for increasing electric field (i) to (iv), with the light polarization in 
the plane of the layer (a) and perpendicular to the plane of the layers (b) 
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7.2.4 Piezo-electric fields in Nitride Quantum wells 

7.2.5 Application: EA modulator 




-1 -2.? -2 -1.J -1 -0..* 
Voltage (V) 

(!) M 

Figure 7.4: High frequency modulation experiments with a 40Gb/s electromodulator. Left: 
absorption versus applied voltage. Right: eye diagram for increasing frequency 




7.2.6 QCSE in Ge quantum wells 

An interesting application of QCSE in a Ge/SiGe quantum well structure. Althought the 
bandgap is indirect, the QCSE effect is strong because of the proximity of the Y point of the 
band structure. 
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Figure 7.5: SiGe sample description 
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Figure 7.6: Band alignement of the various bands of the Ge/SiGe quantum wells structure 
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Figure 7.7: Results: absorption versus applied field, for increasing applied voltages 
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7.3 Cyclotron resonances 
7.3.1 Classical approach 

Let us describe first a purely classical model, where a carrier with charge q and effective 
mass m is subject to a damping term, an a.c. electric field with angular frequency uj and a 

— * 

static magnetic field B along the z-direction: 

dmv mv _ „ „. 

— - + = q(S + vxB). (7.3.8) 

The cyclotron frequency u c is defined as (q has the sign of the particle charge) 



The d.c. Drude conductivity is: 



u c = q — (7.3.9) 

m 



2 

(j = (7.3.10) 

m 



— * 

A tensor of conductivity can be derived from the relationship between j = nqv and the 
electric field (j = a£): 

<~>xx — <~>yy — <J0 r „ 



■ia;)M 
-i 1 



Oxy = -*yx = ^ j-j - ^_ - - (7.3.11) 



Let us consider a now a light with circular (right) polarization; the electric field £ x = —i£ y . 
The conductivity a + is: 



^+ = p" = + ^tt = ff or m _ . = — r 7.3.12 

4 (r m 1 + zw) 2 + ^ t^ 1 + z(u; + u; c ) 

Similarely, the conductivity for the left polarized light writes: 

a = _ 1 ^ r (7.3.13) 

We note that the conductivity cr+ will have a resonance for a; = |c<j c | while the same resonance 
will happen for holes and cr_. The sharpness of the resonance will depend on the product 

UT m . 

7.3.2 Quantum model 

For the quantum model, one starts with the usual Hamiltonian with the substitution for the 
momentum: 

p^(p-qA) (7.3.14) 
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The magnetic field is along the z direction, with a component in the y direction. We choose 
the gauge such that the potential vector is written as: 

A = (Bz sin 9,Bxcos9,0) (7.3.15) 

The total Hamiltonian for a quantum well is then 

H = ^ + V conf (z) + ^(~ih^ + qB ± z) 2 + + q B {l x) 2 + g*» B a ■ B (7.3.16) 

where the component of the field parallel to the growth axis 

B\\ = B cos 9 (7.3.17) 

and the component in the plane 

B ± = Bsin9 (7.3.18) 

have been used. The spin operator has the eigenvalues ±| and the coupling with the magnetic 
field is given by the Lande factor g*. In general, the Hamiltonian 7.3.16 is not solvable 
exactly, and in particular is not separable as the in-plane component of the field is mixed 
with the confinement. In the case that the in-plane component can be neglected, the resulting 
Hamiltonian is: 

H ± = [^ + ^(Hk y + qB ll x) 2 ]<t> n (x) = E n <f> n (x) (7.3.19) 

which, as is apparent, can be mapped onto a Harmonic oscillator. The solutions are the 
celebrated Landau levels with energy: 

E n = h\u c \ (n + -) (7.3.20) 

7.3.3 Determination of the effective mass 

7.3.4 Application: InSb THz detector 

7.3.5 Application: p-Ge laser 

A p-Ge laser is a unipolar laser based on the motion of electrons in crossed electric and 
magnetic fields 



7.4 Interband recombinaison in Magnetic Fields 

In the first order, the bands are shifted by the zero point energy of the Landau levels. The 
situation is much more interesting in confined structures, especially quantum dots where the 
shift of the levels with magnetic field can be used to probe the eigenstates 



Chapter 8 



Second-order processes and 
quasi-particles 

8.1 Interband absorption in indirect materials 

In Si, Ge, AlAs, GaP, the minimum of the conduction band does not occur at the T point, 
while the maximum of the valence band does occur at this location. As a result, interband 
transitions at the bandgap involve always a momentum transfer. This one can be provided 
by any elastic or quasi-elastic process in principle. In practice, however, the most efficient 
process is the emission and absorption of zone-edge optical or acoustical phonons. 




k o k 

Semiconducteur direct Semiconducteur indirect 



Figure 8.1: Schematic band structure diagram of a direct and indirect semiconductor 

The rate at which such process occur can be computed using a second-order, time-dependent 
perturbation theory well described by Feynman diagrams. The rate of such process can be 
written as: 

Rina = \ V I Y^ me j m ^\ H{E c ik l ) - E v (k v ) -hu± E p ). (8.1.1) 

In the above equation, H eo is the coupling between the electron and the incident optical 
radiation; H ep the coupling between the electron and the phonon population. The sum runs 
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over all possible intermediate states The latter is either in the conduction band (if the 
virtual absorption occurs before the emission/absorption of phonon) or in the valence band 
(in the opposite case). The energy 



is conserved in the process, q being the momentum of the phonon. In the latter equation 

1. The 5 function insures that at the end of the process, the energy is conserved. 

2. Two classes of process exist: one in which a phonon is absorbed, one in which one is 
emitted. Obviously only the latter subsists at T = 0. 

3. As the conduction band minimum is at the edge of the Brillouin zone, the momentum 
of the phonon q is much larger than the thermal wavevector of either initial or final 
state. As a result, the process does not conserve the wavevector of the electron or hole. 

4. Therefore, we have to run the integral on both initial an final states with each a 
density of state with an energy dependence Vq(E) = \J — E v and a final state V^(E) = 
■sj 1 E c — E ig . We expect then a energy dependence as (E — E ig ± E p ) 2 . 

Taking also into account the dependence of the matrix elements in the phonon numbers: 



E c (ki) - E v (k v ) -?ku±E p 




and total momentum 



k v + q — k c 




itK 



P va \ 2 \H ac \ 2 N c N v n ph x 2 6(x) + (n ph + l)y 2 6(y) 



(8.1.4) 



8co 



where the phonon occupation number is 



1 




n ph 



exp(^) " 1 



As a result the absorption edge is quadratic energy dependence, and is strongly temperature 
dependence. 
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Figure 8.2: Plot of the square root of the absorption coefficient of Silicon as a function of 
photon energy for various temperatures, as indicated 

The log plot of the absorption of the Germanium shows both the fundamental indirect gap 
at an energy of 0.65-0. 75eV and the fundamental one at 0.85eV. 
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Figure 8.3: Log plot of the absorption of Germanium showing both the indirect fundamental 
gap and the higher energy direct gap 
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Interestingly, the very low absorption (and emissivity) of the indirect materials do not forbid 
a priori the fabrication of relatively efficient emitters. Of course, those can be realized only if 
the minority carrier lifetime can be brought by the reduction of other scattering mechanisms 
(impurities, deep trap, surface states) to a level comparable to the radiative lifetime, equal 
to about 14ms for Silicon at room temperature. Recent results along this line have shown 
efficiencies of more than 1%. However, a laser is believed to be impossible as the cross section 
for absorption is larger than the one for emission. 



8.2 Free-carrier absorption 

Free carriers are responsible for non-resonant light absorption. These phenomena are impor- 
tant parasitic mechanism in optoelectronic devices. 



8.2.1 Classical model 

We had derived in the third chapter the classical response of an electron gas. Assuming this 
classical model with no restoring force 




(8.2.6) 



The imaginary part of the susceptibility is then: 





It can be related to the absorption coefficient a(cu): 



(8.2.8) 
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Figure 8.4: Room temperature free carrier absorption in silicon, for various carrier densities. 

#1: 1.4 x 10 16 cm" 3 #2 8.0 x 10 16 cm~ 3 . #3 1.7 x KFcm" 3 #4 3.2 x KFcm -3 #5 6.1 x 
10 is cm -3 #6 LQ x 10 i9 cm -3 

This relationship predicts a u~ 2 or A 2 dependence of the free carrier absorption; the latter is 
approximately observed in experiments, as shown for the example of Silicon (if one neglects 
an impurity band present in the A = 2 — 5/im )in Fig. 8.4. In other materials, a dependence 
in A p with 2 < p < 3 is observed and is attributed to the variation of the damping term with 
frequency. Note that the specific choice made of the damping term in equ. 3.1.1 implies that 
the Drude scattering time is now r rn = 

As for the case of the interband absorption in indirect materials, the absorption can be 
computed more rigorously using second-order perturbation theory. 

eg yi \ - n \ y (f\H op \n)(n\H ph \i) {f\H ph \n){n\H op \i) 2 

n 

where the sum should be in principle carried over all processes: impurity, optical and acous- 
tical phonons. 



8.3 Excitons 



When a electron-hole pair is created by absorption of a photon inside a semiconductor, the 
pair experiences an attractive Coulomb force. The latter is responsible for the creation of 
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a bound state called an exciton. The exciton can be loosely described by a Hydrogen-like 
system where electrons and holes orbit around their center-of-mass. One distinguish in the 
literature two type of excitons, depending essentially on the ratio of the exciton radius to the 
lattice spacing. When the attraction is so strong that the radius of the exciton is of the same 
order as the lattice spacing, one refers to a Frenckel exciton. Because of the small effective 
mass and the dielectric constant of the semiconductors, the exciton in a semiconductor is 
rather of the Wannier type, where the exciton radius is of many lattice periods. 



8.3.1 Elementary treatement 

Applying an effective mass description of the electron and hole, one easily maps the system 
onto a Hydrogenoid particule with reduced mass fi' 1 = (m*) -1 + (m£) -1 . The energy 
spectrum is the one of the Hydrogen atom: 

E n = -^j (8.3.10) 

rr 

where the renormalized Rydberg energy TZ* is given by 

n* = -^^n (8.3.11) 

m e 2 

were e is the semiconductor dielectric constant. As the energy E n is negative, the exciton 
features will appear as discrete eatures inside the material's bandgap.Both effective mass 
and dielectric constant are strongly dependent on the material's effective mass. As a result, 
and as shown in Fig. 8.5, the binding energy of the exciton grows strongly with the material 
bandgap. In large bandgap materials such as GaN, excitonic features are observed up to 
room temperature, whereas in narrow gap materials such as InAs they are difficult to observe 
even at cryogenic temperature. 
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Figure 8.5: 
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8.3.2 Excitons in confined structures 

Because they are states that cover a large number of lattice spacing, the excitons are well 
described by an envelope function formalism, where the two particule wavefunction is written 
as: 

^(r e ,r h ) = C(k e ,k h )^ ke (r e )tP kh (r h ) (8.3.12) 

However, because the exciton is a localized state, it is more conveniently written in terms of 
a Wannier functions a n (r; Rj) of site Ri rather than in terms of the Bloch functions. Wannier 
functions and Bloch functions are defined by the set of reciprocal relations: 

a n (r; Ri) = — L V] exp(-ikRi)ip nk (r) (8.3.13) 
^ N k 

^nk(r) = —= y^exp(ikRi)a n (r,Ri). (8.3.14) 

Hi 

The exciton wavefunction then writes: 

V(r e ,r h ) = 4r?y2 $ (^> R h )a(r e ; R e )a(r h ; R h ) (8.3.15) 

The exciton wavefunction is solution of the equation 

h 2 h 2 I 

V 2 Te - - — V* - fracq 2 Ane\r e - r h \ 0(r e , r h ) = -E<f)(r e , r h ). (8.3.16) 



2m e e 2m e 

This is a hydrogen-like Hamiltonian that can be solved using the center-of-mass coordonates: 

„ m e r e + m h r h 

R=^^ — (8.3.17) 

m e + m h 

and a relative coordinate 

r = r e -r h . (8.3.18) 

The two resulting equations are 

h * :V 2 R *(R) = E R *(R) (8.3.19) 



2M 

fc2 2 

("2^ V ^ (r) ~ 4^) 0(r) = Er0(r) (8 - 3 - 20) 

where the reduced mass \i and the total mass M = m e + mh have been used. The solutions 
of the center of mass equation are readily obtained: 

V k (r) = -j=exp(-KR) (8.3.21) 

h 2 K 2 

E R = r ±^- (8.3.22) 

2M y 1 
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The solution for the relative motion follows the usual Hydrogen solution, where the wave- 
function is written as a product of a radial function and the spherical harmonics: 

(r) = R n , l (r)Y lim (6,<f>) (8.3.23) 

where R n j are teh Laguerre polynomials and Yi ttn (9, 0) the spherical harmonics. The energy 
spectrum is then given by the same expression as shown above: 

EJn) = E r (oo) - \ (8.3.24) 
n 2 

where the renormalized Rydberg energy is given by: 

TV = «L (8.3.25) 
which the result already obtained in the preceding paragraph. 



8.3.3 Quantum well 

A interesting case is the one of the exciton confined in dimension in a quantum well. The 
envelope function for the exciton 0(r) follows now the Schr odinger equation: 

2 \ / 2 fc 2 7i-^"2 

2^(r) + K//(r)0(r) = (E - — - e[ c) - - E G )J>{r) (8.3.26) 
where the effective interaction potential now writes 

Veff(r) = / d Zl dz 2 - = w — — (8.3.27) 

in the limit of the perfect two-dimensional case, valid when the exciton radius is much larger 
than the quantum well width. The equation is the one of a two-dimentional Hyrdrogen atom. 
The energies of the latter are defined by: 

1Z* 

E n = - (n _ 1/2)2 - ( 8 - 3 - 28 ) 

As a result, the excitonic confinement of the ground state is multiplied by a factor of 4 
compared to the bulk case. As a result, excitonic features are now observable in quantum 
wells at room temperature. 



Chapter 9 
Quantum dots 



9.1 Basic ideas 

We have seen that the gain in semiconductor laser is proportional to the product of the 
density of state and the Fermi distribution difference: 

g{E) = C^(M(/c(M - f v {hw)) (9.1.1) 

where C is a prefactor that contains the matrix element and fundamental constants. Chang- 
ing the dimensionality of the active region, i.e. going from a bulk material to a quantum well 
and finally a quantum wire and dot, because of the change in the density of state, has a 
strong influence on the maximum gain. Proposed already by Arakawa and Sakaki in their 
celebrated paper in 1982, they showed that in general the reduction of dimensionality had 
a very strong influence on lowering the threshold current density and decreasing the tem- 
perature coefficient of the gain. As shown schematically in Fig. 9.1, the density of state is 
narrowed by the reduction of the dimensionality. 




Figure 9.1: Comparison of the density of state for bulk, quantum wells, quantum wires and 
quantum dots.(Asada, IEEE JQE 1986) 
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As shown in Fig. 9.2 the gain should be greatly improved for a given carrier density when 
decreasing the dimensionality of the active region. 



ioV 1 1 1 r 




Wavelength X(jim) 



Figure 9.2: Comparison of the computed gain for bulk, quantum well, quantum wire and 
quantum dot material. (Asada, IEEE JQE 1986) 



The temperature dependence of the threshold should be greatly improved, with a gain inde- 
pendent of the temperature in the limit of the quantum dots, as shown in Fig. 9.3 
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Figure 9.3: Comparison of the temperature dependence of the threshold current for bulk (a), 
quantum well (b), quantum wire (c) and quantum dot material (d).(Arakawa and Sakaki, 
APL 1982) 



9.2 Fabrication issues 



The problem of fabricating quantum dots with good size uniformity, crystalline quality and 
purity is a very difficult technical challenge. A large number of techniques have been used. 
At this point, the most successful is the self-assembly of dots by Stransky-Krastanow mode 
of 3D growth. This technique enables the fabrication of a relatively large density of dots 
with excellent purity and crystallographic quality and with reasonable size uniformity. 




Figure 9.4: Atomic force microscopy image and transmission electron microscopy image of 
self-assembled quantum dots (data courtesy of Andrea Fiore) 
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9.3 Quantum dot luminescence 

The luminescence of self-assembled quantum dots is shown in Fig. 9.5 as a function of pump 
intensity. 




Figure 9.5: Luminescence of self-assembled quantum dots as a function of (data courtesy of 
Andrea Fiore) 

The data show clearly the apparition of successive excited states as the pump intensity 
is increased. Those can be computed by an envelope function computation, taking into 
account strain and the three dimensional confinement. The broadening of the curve is 
almost completely caused by the size variation of the dots. The system is thus completely 
inhomogenously broadened. 

9.4 Quantum dot lasers 

The light versus current characteristics of a laser based on an active region with quantum 
dots is shown in Fig. 9.6. As compared to a quantum well material, the quantum dot has a 
very diffent behavior. The three essential physical differences are: 

1. The active region is inhomogenously broadened 

2. The capture time in the dot and interlevel scattering in each dot is significantly longer 
(about xlO) than in quantum well structures 

3. The total density of state per dot layer is much smaller than the one of a single quantum 
well 

As a result, the devices exhibit in general a much lower gain (but can also show very low 
transparency current for the same reason). They also exhibit spectral hole burning features. 
Finally, the longer capture time and interlevel scattering time is detrimental to the slope 
efficiency and maximum modulation frequency. 



9.4. QUANTUM DOT LASERS 
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Figure 9.6: Optical power versus injected current characteristics of a quantum dot laser, 
(data courtesy of Andrea Fiore) 

However, for the same reasons, these active region may have very interesting applications 
for optical amplifiers and superluminescent diodes. They can also exhibit very low values of 
the linewidth enhancement factor. 



